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SUMMARY 


Upper  and  lower  bounds  on  the  effective  elastic  properties  of 
blaxlally  fiber  reinforced  material  are  derived.  Two  geometric  models 
for  the  composite  have  been  considered,  and  results  are  presented  for 
several  practical  materials.  In  addition,  bounds  on  the  components  of 
the  shell  stiffness  tensor  are  obtained  in  terms  of  the  bounds  of  the 
elastic  moduli. 
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1.  INTRODUCTION 


Theoretical  study  of  the  elastic  behavior  of  fiber  reinforced 
materials  has  to  date  been  devoted  primarily  to  uniaxially  reinforced 
materials.  Hashin  and  Rosen  [1]*  Introduced  the  composite  cylinder 
assemblage  model  as  an  idealization  of  a  uniaxially  reinforced  material. 

By  use  of  variational  methods,  four  of  the  five  independent  elastic  moduli 
were  exactly  calculated  on  the  basis  of  this  model,  and  lower  and  upper 
bounds  were  given  for  a  fifth  modulus.  In  many  cases,  experimental  re¬ 
sults  are  quite  close  to  the  theoretical  predictions  [2,3]. 

The  problem  has  also  been  treated  from  a  very  general  point  of 
view  by  Hill  [4]  and  Hashin  [5].  In  these  Investigations,  no  special  in¬ 
ternal  geometry  is  assumed}  the  fiber  cross  sections  may  have  arbitrary 
shapes  and  their  positions  are  random.  It  should  be  realized  that  under 
these  circumstances  much  less  is  knovm  about  the  internal  geometry  of 
the  composite  than  in  the  case  of  the  composite  cylinder  assemblage  used 
in  [1],  for  in  the  latter  case  the  fiber  cross  sections  must  be  circular 
and  their  arrangement  is  of  a  very  special  nature.  Consequently,  the 
investigations  [4,5]  did  not  permit  exact  calculation  of  any  effective 
modulus,  since  the  general  model  used  did  not  contain  enough  information 
for  unique  determination  of  effective  moduli.  Therefore,  the  approach 
consisted  of  bounding  of  the  effective  moduli. 

The  Investigations  [4,5]  should  be  considered  as  a  first  step 
towards  the  solution  of  a  general  and  difficult  problem;  namely,  estima¬ 
tion  of  effective  moduli  for  Internal  geometry  which  is  specified  in  a 
statistical  manner.  Since  the  bounds  are  based  only  on  very  rudimentary 
geometrical  Information,  i.e,  phase  volume  fractions,  it  is  not  surpris¬ 
ing  that  for  the  elevated  fiber  to  matrix  stiffness  ratios  usually  en¬ 
countered  in  practice,  not  all  the  bounds  give  sufficiently  close  estimates 
to  be  of  practical  importance.  Therefore,  in  what  follows,  the  results 
obtained  in  [1]  will  be  used  for  effective  moduli  of  uniaxially  reinforced 
materials. 


A  number  of  other  investigations  of  the  problem  are  to  be  found 
in  the  literature.  However,  these  are  based  on  very  crude  approximations 
and  thus  do  not  seem  to  be  sufficiently  reliable.  For  a  general  compre¬ 
hensive  survey  of  work  on  properties  of  uniaxially  reinforced  materials, 
see  [6]. 


Practical  experience  and  theoretical  study  of  uniaxially  rein¬ 
forced  materials  reveal  that  these  materials  are  very  stiff  and  strong 
in  fiber  directions.  However,  in  directions  transverse  to  the  fibers, 


*Numbers  in  brackets  refer  to  appended  references. 


both  stiffness  and  strength  are  an  order  of  magnitude  smaller.  In  miny 
applications,  such  as  plate  and  shell  constructions,  stiffness  and 
strength  of  comparable  magnitudes  are  required  in  all  directions.  T'lis 
may  be  achieved  by  constructing  a  material  with  two  or  more  directions 
of  reinforcement.  The  present  investigation  is  concerned  with  the 
theoretical  prediction  of  the  elastic  properties  of  a  biaxially  fiber 
reinforced  material.  The  method  used  is  based  on  bounding  techniques  on 
the  geometrical  model  of  a  fiber  reinforced  material  which  has  been  in¬ 
troduced  in  [1].  It  will  be  seen  that  the  bounds  obtained  are  mostly 
very  close  together  and  thus  give  very  good  estimates  of  the  effective 
elastic  properties. 

The  results  are  then  exploited  for  analysis  of  elastic  behavior 
of  biaxially  reinforced  shells. 
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2.  DEFINITION  OF  EFFECTIVE  ELASTIC  MODULI 


The  material  under  consideration  is  a  biaxially  fiber  rein¬ 
forced  material.  It  is  composed  of  layers  which  are  assumed  to  be  of 
equal  thickness  h.  In  each  layer  the  fibers  are  in  one  direction.  The 
reinforcement  direction  alternates  from  layer  to  layer.  The  material  Is 
referred  to  a  coordinate  system  X2,  x^,  in  which  X2  is  normal  to  the 

layers  and  x^^,  x^  bisect  the  angles  formed  by  alternating  reinforcing 

fibers  (Figure  1).  The  angles  between  fiber  directiot\s  and  the  x^^  axis 

are  (j)  and  (180  -  ({i),  respectively.  It  is  assumed  that  In  each  layer  the 
number  of  fibers  through  the  thickness  is  very  large  and  also  that  in  the 
material  specimen  the  number  of  layers  is  very  large.  The  matrix  and 
fiber  materials  are  assumed  to  be  linear  elastic  and  isotropic.  The 
problem  then  is  to  predict  the  elastic  properties  of  the  biaxially  rein¬ 
forced  m.aterlal  in  terms  of  elastic  moduli  of  the  constituents,  their 
volume  fractions  or  additional  geometrical  characteristics,  and  the  angle 
of  reinforcement  4). 

The  effective  elastic  moduli  of  a  statistically  homogeneous 
elastic  composite  material  are  defined  in  terms  of  stress  and  strain 
averages.  Let  a  large  specimen  be  subjected  to  surface  displacements 
of  the  form 


u°(S) 


(1) 


where  are  constant. 


Then  the  volume  average  strains  defined  by 


■ij 


ij 


are  given  by 


7  «  , 

ij  ij 


(3) 


The  average  stresses  are  necessarily  linearly  related  to  the 
average  strains.  Write 


°ij  “  ^Ijkl  ^kl  “  ^Ijkl 

The  coefficients  are  defined  as  the  effective  elastic  moduli. 

Dually,  let  the  specimen  be  subjected  to  boundary  tractions  of  the  form 

(5) 


T°(S) 


°ij"j’ 


3 


I 


Figure  1.  Geometrical  Model  of  the  Material 
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where  are  the  components  of  the  outward  normal  and  are  constant 
stresses.  Then 


o 


Ij 


(6) 


The  necessarily  linear  relation  between  the  average  stresses  and  strains 
is  now  written  in  the  form 


■ij 


C*  rf° 

ijkl  kl' 


(7) 


The  coefficients  S,,, ,  are  defined  as  the  effective  elastic 

*  * 

compliances.  For  a  very  large  specimen,  the  tensors  and 

are  reciprocal,  and  this  will  be  assumed  forthwith. 


The  effective  elastic  moduli  and  compliances  can  also  be  de¬ 
fined  in  terms  of  stored  elastic  energy.  When  (1)  is  prescribed,  it  may 

be  shown  [1]  that  the  strain  energy  stored  in  a  large  specimen  is 


W 


G 


(8) 


Here  the  specimen  volume  has  been  taken  as  unit  volume  without  loss  of 
generality.  Dually,  when  (5)  is  prescribed,  the  stress  energy  (here  and 
in  what  follows,  elastic  energy  in  terms  of  strains  is  termed  strain 
energy;  elastic  energy  in  terms  of  stresses,  stress  energy)  is  given  by 
[1] 


W 


a 


i.  c*  0  0 

2  ^Ijkl  '’ij  '’kl* 


(9) 


Expressions  (8)  and  (9)  may  be  regarded  as  the  strain  and  stress  energy 
densities  of  the  composites. 


In  general,  the  number  of  effective  elastic  moduli  is  at  most 
21.  Material  symmetry  will  reduce  this  number.  It  is  clear  from  Figure  1 
that  the  present  material  has  three  planes  of  elastic  symmetry,  i.e.,  the 


X 


2  ,  X2  ^3  planes.  Consequently,  the  material  is  ortho¬ 


tropic  and  has  nine  independent  effective  elastic  moduli  or  compliances. 
The  stress-strain  relation  (9)  then  assumes  the  following  form  [7] : 


^11  ^11  ^12  "22  ^13 


(a) 


i 

* 

} 

i 


1 


I 
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N 

CM 

CM 

^12  ^11  S2  *^22  ^23  ^33’ 

(b) 

33  " 

*  _  *  _  *  _ 

^13  ^11  ^3  ^^22  S3  ^33  ’ 

(c) 

(10) 

'12  “ 

^S4  ^12’ 

(a) 

'23  * 

^S5  ^23’ 

(b) 

'31  " 

^Sb  Si’ 

(c) 

(11) 

where  six-by-six  matrix  notation  has  been  used  for  the  It  is 

evident  that  the  coefficients  in  (11  a,b,c)  are  shear  moduli.  Therefore, 
these  will  be  given  the  following  notation: 


44 

* 

55 

* 

66 


'll 


'22 


'33 


^11  '’ll 


^12  '’ll 


^13  '’ll 


* 

“  S2’ 

(a) 

* 

“  S3’ 

(b) 

* 

■Si- 

(c) 

be  similarly  i^itten: 

*  - 

* 

+  Sj^2  ^22 

S3 

Q 

*  _ 

* 

+  S22  O22  + 

S3 

Ojy 

*  _ 

* 

+  S23  022 

^33 

'’33, 

(a) 

(b) 

(c) 


(12) 


(13) 


'12 


'23 


'31 


2S44 


2S55  023, 


^^66  '’31* 


(a) 

(b) 

(c)  (14) 
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*  * 

The  three-by-three  S3nninetrlc  matrices  and  S^j  (1,  j 


1,  2,  3)  are  reciprocal. 

(13)  In  terms  of  Young's  moduli  and  Poisson's  ratios  [7] 


It  Is  common  practice  to  write  the  In 


Thus, 


11 


1. 

* 


(a) 


22 


(b) 


33 


*  » 

"3 


(c) 


(15) 
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12 


23 


31 


V 


E 


'1 

* 

23 

* 

2 

* 

31 

* 


* 

7 

'^2 

* 

* 

h 

* 

* 

? 

"1 


.  (a) 


=  -  ^  ,  (b) 


»  .  (c) 


(16) 


To  explain  the  meaning  of  these,  assume,  for  example,  uniaxial 
average  stress  and  all  others  vanish.  Then 


^  -  -  ^13  - 

^22  *  "  *  °11’  ^33  “  "  *  '^11’ 


and  the  general  stress-strain  law  (13)  can  be  obtained  by  superposition 
of  the  normal  stresses  in  the  coordinate  directions. 
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3.  GEOMETRICAL  MODELS  OF  BIAXIALLY  FIBER  REINFORCED  MATERIAL 


For  the  purpose  of  prediction  of  effective  elastic  properties, 
two  geometrical  models  of  the  composite  will  be  considered.  The  first  of 
these  is  termed  the  layer  model  (LM) ;  and  the  second,  the  composite 
cylinder  model  (CCM) . 

LAYER  MODEL 


The  material  is  composed  of  homogeneous  elastic  layers  of  equal 
thickness.  Each  layer  has  the  elastic  properties  of  a  uniaxially  fiber 
reinforced  material.  Such  materials  are  transversely  isotropic,  the 
fiber  direction  being  the  axis  of  elastic  symmetry.  As  a  result  of  this 
symmetry,  such  a  material  has  five  independent,  effective,  elastic  moduli. 

It  should  be  borne  in  mind  that  in  order  to  bring  out  this 
symmetry,  it  is  necessary  to  use  a  coordinate  system,  one  of  whose  axes 
is  in  the  fiber  direction  while  the  other  two  are  normal  to  it.  These 
coordinate  systems  are  shown  in  Figure  2,  with  respect  to  the  two  differ¬ 
ent  fiber  layers. 

In  terms  of  the  model  described  above,  the  problem  takes  on  the 
following  form:  given  the  effective  elastic  moduli  of  a  layer,  or  bounds 
thereon,  with  respect  to  the  primed  and  double-primed  coordinate  systems, 
find  the  effective  elastic  moduli,  or  bounds  thereon,  for  the  layer 
composite. 

COMPOSITE  CYLINDER  MODEL 


In  this  model,  more  detailed  account  is  taken  of  the  internal 
geometry  of  the  composite  than  in  the  LM.  It  is  assumed  that  each  layer 
is  described  by  the  composite  cylinder  assemblage  which  has  been  intro¬ 
duced  in  [1].  To  describe  such  an  assemblage,  assume  that  one  has  at 
his  disposal  composite  cylinders,  each  of  which  consists  of  a  circular 
cylindrical  fiber  and  a  concentric  circular  matrix  shell,  perfectly 
bonded  together.  In  each  cylinder,  the  ratio  between  fiber  radius  and 
outer  matrix  shell  ratio  is  the  same  (Figure  3(a)).  The  volume  of  each 
layer  is  imagined  to  be  filled  out  progressively  with  parallel  composite 
cylinders  of  decreasing  outer  radii,  (Figure  3(b)).  In  the  limit,  the 
whole  volume  is  thus  filled  out,  and  the  composite  analysis  then  becomes 
based  on  the  analysis  of  a  single  composite  cylinder. 
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Figure  2.  Reference  Frames. 


a/b  s  constant 

(a)  COMPOSITE  CYLINDER 


Figure  3.  Composite  Cylinder  Model. 


A.  VARIATIONAL  BOUNDING  METHODS  FOR  EFFECTIVE  ELASTIC  MODULI 


As  has  been  shown  in  Section  2,  the  determination  of  effective 
elastic  moduli  is  equivalent  to  the  determination  of  stress  averages  for 
a  boundary  condition  of  type  (1),  or  of  strain  average  for  a  boundary 
condition  of  type  (5).  In  general,  in  order  to  find  an  average,  the  de¬ 
tailed  field  must  first  be  computed.  It  is  clear  that  field  determina¬ 
tion  in  biaxially  fiber  reinforced  materials  is  an  Intractable  problem. 

This  difficulty  can  be  avoided  to  a  large  extent  by  variational 
bounding  techniques.  Here,  the  aim  is  to  find  lower  and  upper  bounds  for 
the  effective  moduli  and  compliances,  instead  of  attempting  to  compute 
them. 


The  possibility  of  variational  bounding  hinges  upon  the  results 
(8)  and  (9),  and  the  variational  principles  of  the  theory  of  elasticity. 
According  to  (8)  and  (9),  the  strain  or  stress  energy  is  given  in  terms 
of  strains  or  stresses  known  from  the  boundary  conditions,  and  the 
effective  elastic  moduli  or  compliances.  The  elastic  energy  can  be 

bounded  by  use  of  the  variational  principles  of  elasticity,  and  thus 

*  * 


bounds  are  obtained  on  combinations  of  C 


ijkl 


or 


^Ijkl* 


BOUNDING  METHOD  FOR  LAYER  MODEL 


Let  a  displacement  field  of  type  (1)  be  prescribed  on  the 
boundary.  An  admissible  displacement  field.  u^(x),  for  the  principle  of 

minimum  potential  energy  (see  [8],  for  example)  is  defined  as  a  continuous 
field  which  satisfies  the  boundary  conditions.  The  simplest  admissible 
field  for  the  present  case  is 

u^(x)  ■  ^ij^J» 

which  obviously  satisfies  all  the  requirements. 

Now  compute  the  integral 


7  /c 


ijkl  ^Ij 


^kl 


(18) 


throughout  the  LM  specimen.  Here, 


e 


Ij 


(19) 


and 

C,  ,  assume 
ijkl 


are  the  space  variable  elastic  moduli.  In  the  present  case 
only  two  different  values'  because  there  are  only  two  kinds 
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of  layers.  Denote  these  by  for  the  layers  in  which  the  fibers  are 

inclined  at  angle  (j)  with  respect  to  and  by  for  the  layers  in 

which  the  fibers  are  inclined  at  angle  (180  -i(>)  with  respect  to 

(Figure  2).  "*^ijkl  elastic  moduli  of  the  transversely 

isotropic  layers  referred  to  in  the  specimen  coordinate  system  x^^,  X2,  x^. 

Since  the  layer  moduli  are  defined  with  respect  to  the  coordinate  systems 
Xp  Xp  or  Xp  Xp  x^  (Figure  2),  the  values  of  and  are 

defined  by  fourth-rank  tensor  transformation  rules. 

Now  in  the  present  case,  from  (17)  and  (19), 


Assuming  equal  volumes  of  different  layers,  and  substituting  (20)  into 
(18),  the  Integral  (18)  reduces  to 

2W  ~  j  f'^ijkl  "^ijkl^  ^ij  ^kl  • 

From  the  principle  of  minimum  potential  energy, 

.  (22) 

Inserting  (21)  and  (8),  respectively,  into  (22),  one  finds 

^ijkl  ^ij  ^kl  -2  ^'^ijkl  "^ijkl^  ^ij  '^kl  * 

One  can  now  proceed  in  like  manner  with  the  specimen  loaded  by 
tractions  of  type  (5) .  Here  an  admissible  stress  system  is  needed  in 
conjunction  with  the  principle  of  minimum  complementary  energy  [8].  An 
admissible  stress  system,  d^j(2i)»  is  defined  as  one  which  satisfies 

equilibrium,  traction  continuity,  and  the  boundary  conditions.  The 
simplest  9^^  in  the  present  case  is  the  constant  stress  system. 


Now  compute  the  Integral 


where  S^jj^j^(2c^)  are  the  space  variable  compliances,  which  again  assume 
only  two  values  and  in  the  different  layers.  They  are  the 


reciprocals  of  J^®spectlvely .  Proceeding  as  before 

for  unit  volume,  (25)  assumes  the  form 


2W' 


a  M  1  <'  • 


2  ^'^IJkl  "^Ijkl^  °ij  '^kl  ‘ 


Because  of  the  principle  of  minimum  complementary  energy. 

Inserting  (9)  and  (26)  into  (27),  it  is  found  that 

* 


o  o  0 

^IJkl  °1J  ‘^kl 


i  2  "\jkl>  “IJ  “kl 


(26) 


(27) 


(28) 


BOUNDING  METHOD  FOR  COMPOSITE  CYLINDER  MODEL 


Let  the  volume  filled  out  by  composite  cylinders  in  the  CCM  be 
denoted  by  V_  and  the  remaining  volume  by  Vj^.  Thus,  by  hypothesis, 
may  by  made  indefinitely  small  and  in  the  limit 

V.  (29) 

Let  the  biaxially  reinforced  specimen  be  subjected  to  a  dis¬ 
placement  of  type  (1)  on  the  boundary.  An  admissible  displacement  field 
is  now  conctructed  as  follows:  The  field  (17)  is  Imposed  throughout 
Vj^.  The  composite  cylinders  are  then  subjected  to  (17)  on  their  bound¬ 
aries.  This  defines  an  elasticity  boundary  value  problem  for  a  composite 
cylinder.  Since  by  construction  all  composite  cylinders  are  geometrically 
similar,  the  same  solution  is  valid  for  each  of  them.  Let  the  displace¬ 
ment  in  a  composite  cylinder  thus  obtained  be  denoted  by  u^j^.  Then  the 
admissible  displacement  field  for  the  composite  is 

-  u^  in  V^,  (a) 

“l  "  ^ij  *j 

Obviously  this  field  satisfies  continuity  and  (1)  on  the  boundary  and  is 
thus  admissible. 

Now  (18)  is  computed  using  (30).  The  result  can  be  written  in 

the  form 


W‘ 


(31) 
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where  W^n)c  strain  energy  in  the  nth  composite  cylinder  subjected 

to  (1)  on  its  boundary  and  is  (18)  for  V|^.  Since  the  strains 
associated  with  (b)  are  constant  is  proportional  to  Vj^.  Since 

Vf^  disappears  in’  the  limit,  so  does  *W^r*  Thus,  in  the  limit 


W 


e 


I 

n  “  1 


(32) 


Because  of  the  geometrical  similarity,  the  strain  energy 
'^%)c  composite  cylinder  can  be  written  in  the  form 


W 


(n)c 


=  w 


V 


(n)c 


(33) 


where  w^  is  one  and  the  same  number  for  all  cylinders  and  is  a  function 
of  ,  the  elastic  moduli  of  fibers  and  matrix  and  the  ratio  between 
the  radii  a/b. 


Therefore,  (32)  assumes  the  form 


n  -  N  -+■  “ 

-  w^  [  V 
n  ■  1 


(n)c 


w^V. 


(34) 


Now  from  (8) ,  (34) ,  and  the  principle  of  minimum  potential  energy, 


1 

2 


"kl  i.'' 


(35) 


where  V  has  been  taken  as  unit  value.  It  should  be  remembered  that  in 
order  to  find  w^,  only  one  composite  cylinder  has  to  be  analyzed. 


A  similar  procedure  is  now  used  with  the  admissible  stresses. 
For  boundary  conditions  of  type  (5),  the  stress  system  is  Imposed 
on  Vj^.  Then  the  composite  cylinders  are  subject  to  the  tractions 


^i  “  '"ij 

(36) 

on  their  boundaries.  The  resulting  stresses 
ders  are  •  The  admissible  stress  system 
is  then 

Inside  the  composite  cylln- 
In  the  composite  specimen 

Oij  -  ojj 

in  V^, 

(a) 

°ij  “  '^Ij 

in  V^. 

(b)  (37) 
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Computation  of  (25)  with  (37)  yields 


W*"  -  I  W?„)c  +  W°R  , 


(38) 


where  ^n)c  strain  energy  in  the  nth  composite  cylinder  sub¬ 

jected  to  (36)  on  its  boundary.  Again  in  the  limit  as  Vj^  becomes  in¬ 
finitesimal,  so  does  .  Then  by  reasoning  completely  analogous  to 
that  given  above, 

-  W°  V,  (39) 

where  is  defined  by 


-  w<^  V 


(n)c‘ 


(40) 


Application  of  minimum  complementary  energy  with  (39)  and  (9), 
for  unit  composite  volume,  yields 


2  Sijkl  Hj 


4l 


(41) 


It  turns  out  that  the  variational  bounding  methods  which  have 
been  described  can  bg  advantageously  exploited  only  for  bounding  of  the 
shear  moduli  G*2  ,  G23  ,  G31  ,  defined  in  (11)  and  (12).  The  reason 
for  this  is  that  with  the  present  orthotropic  stress-strain  laws  (10), 
(11),  (13),  and  (14),  the  energies  on  the  left  sides  of  (33),  (28),  (35)> 
and  (41)  can  be  expressed  in  terms  of  one  single  effective  modulus  or 
compliance,  only  for  the  case  of  an  average  pure  shear  strain  or  stress. 
This  is  because  only  the  shear  stress-strain  relations  (11)  and  (14)  are 
given  in  terms  of  one  single  modulus  or  compliance.  In  these  circum¬ 
stances,  a  bound  on  energy  gives  a  bound  on  an  effective  modulus  or 
compliance.  However,  in  the  case  of  average  normal  strains  or  stresses, 
the  energies  on  the  left  sides  of  (23),  (28),  (35),  and  (41)  involve 
several  moduli  or  compliances.  Thus,  bounds  on  combinations  of  these  are 
obtained,  and  it  is  a  very  difficult  matter  to  deduce  from  these  bounds 
on  the  individual  moduli  and  compliances.  These  difficulties  will  be 
further  explained  below.  Consequently,  the  variational  bounding  methods 
as  described  here  will  not  be  used  for  estimation  of  normal  effective 
moduli  and  compliances. 
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5.  RESUME  OF  ELASTIC  PROPERTIES  OF  UNIAXIALLY  FIBER  REINFORCED  MATERIALS 

Knowledge  of  the  clastic  properties  of  unlaxlally  fiber  rein-  | 

forced  materials  Is  of  fundamental  Importance  for  the  present  study,  since  | 

the  blaxlally  reinforced  composite  Is  composed  of  unlaxlally  reinforced 
layers.  Indeed,  every  effort  will  be  made  to  base  prediction  of  the 
elastic  moduli  of  the  blaxlally  reinforced  material  on  the  moduli  of  a 
unlaxlally  reinforced  material. 

Unlaxlally  fiber  reinforced  materials  are  transversely  Iso¬ 
tropic,  the  fiber  direction  being  the  axis  of  transverse  Isotropy.  Let 
the  fiber  direction  be  xj^,  while  X2,  X3  are  the  other  two  axes  of  the 
Cartesian  reference  system.  Figure  2a.  Similarly,  an  x^,  X2,  X3  system 
of  axes  may  be  taken  for  the  other  layer.  Figure  2b. 

The  material  has  five  Independent,  effective  elastic  moduli  and 
the  effective  stress-strain  law  may  be  written  in  the  following  form  [1]: 

(a) 

(b) 

(c)  (42) 


(a) 

(b) 

(c)  (43) 

Here,  the  overbar  denotes  average  over  the  layer  volume.  The 
average  strains  may  be  expressed  In  terms  of  average  stresses  by  rela¬ 
tions  similar  to  (42)  and  (43),  via  the  effective  elastic  compliances. 

elastic  moduli  of  technical  Importance  will  now  be 
defined  In  terms  of  the  moduli  In  (42)  and  (43): 


11 


’22 


-  C, 


-  C, 
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11  ^12  ^22 

1  ^  1 
S2  *^33’ 

.1  t 

I 

11  ^22  ^22 

S3  ^33’ 

1  »  _t 

?  _  • 

11  S3  ^^22 

S2  ^33* 

_i  »  _i 

^12  "  ^^4  '^12’ 

_  *  »  » 

I 

S3  ■  ^^22  '  S3^  ^^23’ 

_i  1  _i 

Si  "  ^^4  ^31* 

Kt  -  2  (^22  +  C23), 


L5 


Plane  strain  bulk  modulus 


(44) 


-^0^  «■»> . » 


Transverse  shear  modulus 


Axial  shear  modulus 

Axial  Young's  modulus 


Axial  Poisson's  ratio 


Transverse  Young's  modulus 


where 


Transverse  Poisson's  ratio 


Also,  note  the  relation 


G  -  7  (c'  -  C23),  (45) 


t 


(46) 

(47) 


^  ^t  ^ 

^  “  Kj.  +  Gj.’ 


-  1  + 


4  K  V 
t  a 

E  » 
a 


^  ^  ^t 

't  "  +  ip  Gj.' 


(49  a) 


(49  b) 


(50) 


®t  "  2  (1  +  vj' 


(51) 


which  follows  from  (49)  and  (50)  and  Is  of  the  same  form  as  an  Isotropic 
elasticity  relation. 

The  physical  significance  of  the  moduli  (49)andC50)  will  now 
be  briefly  summarized. 

The  modulus  (44)  corresponds  to  the  plane  strain  situation 


til  -  0,  C22  “  E33, 


^12  "  ^23  "  "31  ’  °* 


(52) 
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(53) 


Then,  from  (42  b,  c)  and  (44) > 

- '  -  *  -  *  - '  * 

°22  °33  "  ^  ^^22  ^^33^' 

The  significance  of  (45)  and  (46)  is  obvious  from  (43  b)  and 
(^^  »  rsspectlyelY*  The  shear  modulus  Gt  is  associated  with  shear 

in  the  transverse  X2  X3  plane,  normal  to  the  fibers.  On  the  other  hand, 
the  shear  modulus  is  associated  with  shear  either  in  the  x'  x'  plane 
or  the  xi  X3  plane.  ^  ^ 

The  Young's  modulus  (47)  corresponds  to  the  situation 

. 

_  f  I  _• 

°22  "  °33  "  ®11  ^ 


“12  ■  “23  ■  “31  ■  °- 


(54) 


which  represents  uniaxial  stress  in  fiber  direction. 
Then. 


°11 

11  "  E  * 

d 


(55) 


Also,  in  the  case  (54), 


*^22  "  ^33  “  "  %  *^11’ 


(56) 


which  defines  the  Poisson's  ratio  (48). 


The  Young's  modulus  (49  a)  is  associated  with  either  one  of 
the  loadings 


or 


o'  ■  oln  ■  0, 

o' 

0  > 

11  22  ’ 

33 

^  0 

j 

■  “l2- 

“Ii  -  “33  ■ 

CM 

CM 

D 

Then  for  (57  a). 

-  °33 

r  * 

c-  1 

*^33  "  ’ 

^22 

-  -v^. 

C33  , 

-  o'. 


“a"  “23  ■  “a  ■  “■ 


(a) 

(b)  (57) 


(a) 
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■an 


WVrt*'—  r- 


I 

K 


-V 


t 


which  also  define  the  Poisson's  ratio  (50), 


(b) 


(58) 


terest. 


In  addition,  the  modulus  in  (42  a)  is  of  some  physical  in- 
It  is  associated  with  the  situation 


r  *  ■  £  * 

^22  33 


0, 


(59) 


which  la  the  case  of  a  specimen  unlaxlally  loaded  while  transverse  de¬ 
formation  is  prevented  by  a  rigid  bonded  enclosure. 


Known  expressions  for  effective  elastic  moduli  of  a  unlaxlally 
fiber  reinforced  material,  analyzed  on  the  basis  of  the  composite  cylin¬ 
der  model  (CCM)  will  now  be  given.  Quantities  for  matrix  material  will 
be  given  the  subscript  m,  and  quantities  for  fibers  will  be  given  the 
subscript  f.  The  matrix  fibers  are  assumed  to  be  elastic  isotropic.  The 
volume  fractions  are  v^,  Vf ,  Vjj  +  V£  -  1,  Isotropic  elastic  moduli  used 
are  C  -  shear  modulus,  v  -  Poisson's  ratio,  e  -  Young's  modulus,  and 
k  *  X  +  G  -  plane  strain  bulk  modulus. 


G 

a 


E 


a 


V 


a 


^  '^m  ^f  ^''m  ”  ''f^ 

E“  Vo  +  Ef  Vj  +  (v  /k  )  +  (v  /k  )  +  (1/G^)  » 

m  I  I  m  m 

V  V,  (v-  -  V  )  (i-  -  f-) 

.  .  m  f  f  m  km  kf 

''o  ^m  +  %  ^m  +  (v^/kf)  +  (Vf/k^)  +  (1/G„)  * 


(60) 

(61) 

(62) 

(63) 


An  excellent  numerical  approximation  to  (62)  is  obtained  by 
omission  of  the  last  term.  Generally,  omisElon  of  the  last  term  in  (63) 
also  yields  a  good  approximation. 
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The  modulus  G^.  has  been  bounded  from  above  and  below.  It  is 

Here  neceesar,.  in  view  ’ho Jali^de^oiefr 

two  cases.  For  one  composite  cylinder,  me 

and  respectively.  For  the  fiber  reinforced  material  re- 

g^ded  ae  a  compoaite  cylinder  aasen.bl.ge,  the  bounda  are  denoted  by 


0<->  and 


Then, 


.(-)  . 


1  +  V, 


.  f  m 


(k  +  2  G^)  V  ' 
^  m _ m  m 

2(k  +  GJ 
m  m 


“°t(c)- 


Equation  (64)  is  a  general  bound  for  G^  for  arbitrary  fiber 
geometry  [5].  Therefore,  it  is  also  a  valid  lower  bound  for  the  CCM. 
However,  alnce  it  can  be  shown  that  G<->>  (64)  la  a  better  lower 

bound  than 

The  expression  for  is  given  by  (138)  in  Section  9. 

With  the  bounds  (64)andC65)  and  expressions  (60),  (62),  and 
(63) ,  follow  bounds  for  and  (49)  and  6o) .  It  should  be  carefu  y 

noted  that,  in  general, 

_(+)  _  g  (G^^^),  ®  (6i 


(G^“h. 


However , 


Here,  the  right  aides  in  (66)  and  (67)  express  the  functional  relations 
defined  by  (49)  and  (50) . 
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6.  TRANSFORMATION  FORMULAE 


It  will  be  recalled  that  the  stress-strain  relations  (10), (11), 
(13),  and  (LA)  of  the  biaxially  reinforced  material  hold  in  the  X2,  x^ 

system  shown  in  Figure  1,  while  the  stress-strain  relations  (A2)and(A3) 

of  the  unlaxially  reinforced  layers  hold  in  the  xj^,  x^,  x^  or  Xp  x'^,  x'^ 

systems  shown  in  Figure  2,  Therefore,  it  is  necessary  to  have  transforma¬ 
tion  formulae  for  stresses,  strains,  and  elastic  moduli  from  one  system 
to  another. 

STRESS  AND  STRAIN  TRANSFORMATION 


Define  the  direction  cosines 


“11 

■ 

cos 

(xj. 

xp 

cos 

4 

• 

“22 

■ 

cos 

(xj, 

X2) 

1  . 

“33 

■ 

cos 

(X3, 

X3) 

■ 

cos 

4 

• 

“12 

m 

cos 

(x{, 

X2) 

■ 

0  . 

“21 

m 

cos 

(xj. 

Xi) 

■ 

0  , 

“23 

m 

cos 

(xj. 

X3) 

• 

0  , 

“32 

m 

cos 

(x^, 

X2) 

■ 

0  . 

“31 

m 

cos 

(X3, 

Xi) 

m 

-  sin 

4  , 

“13 

- 

cos 

(xj, 

X3) 

- 

sin 

4 

• 

(68) 


To  transform  stresses  or  strains  in  the  x^^,  x^,  x^  system 
into  stresses  or  strains  in  the  xj^,  x*,  x^'  system,  use  the  tensor 
transformation  rules 


o'  -  a.  a,  0  , 

ij  ir  j.s  rs’ 

(a) 

IJ  ir  js  rs 

(b) 

(69) 

Substituting  (68)  into  (69),  it  follows  that 


o  *  ■ 

11 


2  2 

^ii  cos  41  +  0^2  sin  4+2  sin  4  cos  4i  (a) 
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Oil  -  cos^  ^  +  0^^  8ln^<^  +  2  sin  ^  cos  <p  , 

033  -  o[[  sln^  ^  +  0^^  cos^  ^  2  °l3  *  cos  , 

2  2 

aj^3  -  -(o[|  -  O3P  sin  41  cos  (|>  +  o^^Ccos  (},  -  sin  (p), 
°22  "  °22  * 

0j^2  “  ■  ^12  “  ^^23 

O23  -  sin  41  -  0*^  cos  (p  , 

with  similar  expressions  for  strain  transformation. 


(a) 

(b) 

(c) 

(d) 

(e) 

(f)  (73) 


ELASTIC  MODULI  TRANSFORMATION 


In  view  of  equations  (23)  and  (28),  it  is  necessary  to  trans¬ 
form  the  elastic^  moduli  of  the  transversely  Isotropic  layers  from  the 
X]^,  X2,  X3  and  Xj^,  X2,  X3  systems  to  the  Xj^,  X2,  X3  system. 

Denote  the  elastic  moduli  In  (42)  and  (43)  by  These  are 

the  moduli  of  a  transversely  Isotropic  layer  with  respect  to  its  own 
(primed)  system  of  coordinates,  which  is  the  only  coordinate  system  which 
brings  forth  the  transverse  symmetry.  Similarly,  the  moduli  of  the  layer 
whose  fibers  are  at  (180  -  41)  orientation  may  be  denoted  as 

It  is  obvious  that 

■  'Ijkl-  (w 

since  the  layers  are  made  of  the  same  material. 


denoted 


Let  the  Cl ,  moduli, 
ijkl 

by  Then, 

'C,,..  -  o', 
ijkl  pi 


transformed  to  the  x^^, 


C 

pqrs 


P 


x^,  X3  system  be 


(75) 


where  momentarily  a!,  are  the  direction  cosines  (68).  Similarly,  trans¬ 
form  Cjjkl  ‘^he  x^i  X2,  X3  system  and  denote  the  result  by  Then 


llr  _  n,  »  *  />.*  *  rf 

u...  1  "o.a.a,  o,C 
ijkl  pi  qj  rk  si  pqrs 


(76) 


Here,  are  the  direction  cosines  (68)  with  ({>  replaced  by 
(180  -  (j))  and  (74)  has  been  used.  The  transformed  moduli  'C,  ,  and 
^ijkl  (23). 
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In  view  of  (42)  and  (43),  the  have  the  following  nonvanishing  com¬ 
ponents:  ^ 


t 

■ 

c' 

1111 

^11  ’ 

1 

■ 

C' 

■ 

c*  . 

2222 

^3333 

22  ’ 

f 

■ 

C' 

a 

c 

■  cl«  * 

1122 

^2211 

1133 

3311 

t 

■ 

c* 

m 

C'  . 

2233 

U3322 

23’ 

t 

■ 

c' 

m 

C' 

“Cl  „  “ 

1212 

2121 

1221 

2112 

f 

• 

C* 

m 

C' 

“  C'  “ 

1313 

^3131 

^1331 

3113 

t 

• 

c* 

m 

C' 

_  r*  ■ 

2323 

^3232 

1^2332 

3223 

^2323  “  ^3232  “  ^2332  “  ^3223  “  2  ^^22  ”  ^23^* 

Considering  now  the  transformation  laws  (75),  (76)  with  (68),  the 
following  features  may  be  observed:  Whenever  subscripts  1,  2,  3,  in 

and  appear  only  even  numbers  of  times,  the  transformation 

laws  (75)  and  (76)  involve  only  even  powers  of  sin  and  cos  ij)  or  of 
sin  (180  -  4i).  Therefore,  in  such  cases,  both  transformation  laws  give 
the  same  results.  Thus,  the  following  relations  result: 

^1111  ^1111  ’  1212  ^1212  » 


Ip  ■  * ' f 

^2323  2323  ’ 


'1122  » 


I P  «  I '  r 

^2233  2233, 

ip  ■  '  T 

^3311  3311’ 


On  the  other  hand,  in  all  the  remaining  and  cos 

and  cos  (180  -  ij))  appear  only  in  odd  powers.  Therefore,  in  this  case. 


and,  therefore,  these  components  do  not  enter  into  (23).  Consequently, 
the  elastic  moduli  tensor 
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(79) 


^ijkl  “2  "*^ijkl) 

has  only  the  components  (78),  Accordingly.  is  orthotropic. 

Obviously,  completely  analogous  reasoning  applies  for  the 
compliances.  Accordingly, 

Sljkl  “  7  ('Sijkl  +  "Sijkl)  (®0) 

Is  orthotropic. 

The  components  of  Cj^jkl  will  now  be  given  explicitly  In  terms  of 
C|jki  on  the  basis  of  (75)  and  (68).  Six-by-six  matrix  notation  will 
again  be  used;  l.e., 

^1111  ”  ^11*  ^1122  ”  etc., 

^1212  "  ^44’  ^2323  “  ^55»  ^3131  “  ^66* 

A  long  and  tedious  calculation  yields 

A  9  0 

■  cos  <p  +  ain^<ji  C22  2  sin  ij)  cos  (P  0^2 


+  4  sin2(;>  cos'‘(;>  (a) 

C  22  ■  C^2  *  (b) 

^22  ■  sin^^i  CJ2  +  co8^(ji  C^2  2  sin^(J)  cos^(()  0^2 

+  4  sin^(p  C03^(p  (c) 

Cj^2  ■  cos^({i  Cj^2  8in^(}!  (d) 

C23  -  sin^^i  0^2  +  cos2(j.  0^3,  (e) 

C{3  -  co82(fr  sin^(;i  +  C08^(j)  8in^(j)  C22 
+  (sin^^i  +  cos^4i)  0^2  “  ^  sln^iji  cos^4>  (f) 

-  cos^^i  C44  +  sin^^i  j  (C22  "  ^23),  (g) 
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(h) 


=  8in2(()  Cl^  +  cos^4>  j  (C^2  “  ^23^* 

Cgg  =  cosset)  sin2(^  (Cjj^  +  CJ2  -  2  C{2) 

+  (sin^(|)  -  cos^(|))^  Cj^  .  (i)  (82) 

It  should  be  remembered  that  the  primed  moduli  on  the  right- 
hand  sides  are  defined  by  (42)  and  (43). 

Similar  results  can  be  written  for  compliances. 
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7.  BOUNDING  FOR  LAYER  MODEL  (LM) 


The  fundamental  relations  for  bounding  of  effective  elastic 
moduli  and  compliances  on  the  basis  of  the  LM  are  (23)  and  (28),  which 
yield 


(C 


ijkl 


(S 


ijkl 


C 

S 


*  X  o  O  - 

ijkl^  ^ij  ^kl  -  °  . 

(83) 

*  \  o  O  f. 

Ijkl^  ‘^ij  °kl  -  °  • 

(84) 

Here,  and  are  known  matrices  which  are  given  in 

terms  of  the  known  single-layer  moduli  by  the  transformation  (82),  where- 
as  and  are  the  unknown  composite  moduli  and  compliance  mat¬ 

rices.  Thus,  (83)  and  (84)  are  equivalent  to  the  statement  that  the 
matrices  and  are  positive  definite. 


It  will  now  be  convenient  to  switch  to  six-by-six  matrix  nota¬ 
tion.  For  moduli  and  compliances,  the  six-by-six  notatioii  Is  defined  by 
(81).  The  strains  and  stresses  are  redefined  according  to  the  scheme 


Ell  - 

^22  “ 

^2* 

633  -  e 

2^2  •  "4* 

2023 

-  S 

’  2^3* 

Oil  “  ®1’ 

®22  " 

'•3o  " 

•k 

■ 

CM 

D 

023  “ 

05. 

^^31  "  ° 

Also,  denote  for  ccnveni  .'f 
^^ijkl  “  ^Ijkl  ■  ^ijkl  » 

^^Ijkl  “  ^Ijkl  “  ^l.ikl  * 


(85) 


(86) 


(a) 

(87) 


It  should  be  remembered  that  all  the  matrices  in  (87)  are 
orthotropic.  Thus,  these  matrices  written  out  in  detail  in  six-by-six 
matrix  notation  are 
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^^11 

^^12 

"^13 

0 

0 

0 

^^2 

^^3 

0 

0 

0 

ACi3 

AC23 

AC33 
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0 

0 

0 

0 

0 

AC,, 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^^66 

^^11 

AS_ 

.1^ 

"^13 
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0 

0 

^'l2 

^^22 

^^23 

0 

0 

0 

AS33 

AS23 

AS33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^^66 

- 

-  T 

(88) 


(89) 


The  matrices  (88)  and  (89)  must  be  positive  definite  according 
to  (83)  and  (8A).  It  is  seen  that  these  matrices  are  completely  defined 
by  the  submatrices. 


AC 


AC 


^^11 

ACi3 

^^12 

^^23 

f 

-13 

"^23 

(90) 

0 

0 

0 

0 

> 

0 

0 

(91) 
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^^11 

^^2 

^^13 

^^22 

"^23 

» 

^13 

AS23 

AS33^ 

(92) 

_ 

^^44 

0 

0 

0 

0 

• 

0 

0 

(93) 

For  AC  to  be  positive  definite,  it  is  necessary  and  sufficient 
1  2 

that  A£  and  AC  be  positive  definite.  Similarly,  for  AS  to  be  positive 

1  2 

definite,  it  is  necessary  and  sufficient  that  AS  and  AS  be  positive 
definite. 

As  is  well  known,  a  necessary  and  sufficient  condition  for  a 
matrix  to  be  positive  definite  is  that  Its  eigenvalues  be  positive.  Since 
2  2 

the  matrices  AC  and  Afi  are  already  in  diagonal  form,  as  can  be  seen 
from  (91)  and  T93) ,  the  required  positive  definiteness  of  these  matrices 


conditions 

^^44  -  ®  » 

(a) 

®  ' 

(b) 

(c) 

(94) 

^^44  >  0  , 

(a) 

^^55  ^  0  » 

(b) 

^^66  >  0  , 

(c) 

(95) 

In  view  of  (87  a),  in  six-by-six  notation,  (94)  can  be  written 
in  the  form 


^44  C4A  * 


^55  -  Ss 

Se  -  ^66 


(a) 

(b) 


(c)  (96) 


which  provides  upper  bounds  for  the  shear  moduli  of  the  biaxially  re¬ 
inforced  material. 

In  view  of  (87  b),  in  six-by-six  matrix  notation,  (95)  assumes 

the  form 


^44  -  ^A4  ’ 
^55  -  ^55  * 


(a) 


(b) 


hethi  ■ 


(c) 


(97) 


Now  recall  that  since  and  are  effective  elastic 

moduli  and  effective  compliances  of  the  biaxially  reinforced  material, 
they  are  reciprocal  to  one  another.  It  should  also  be  remembered  that 
because  of  orthotropy,  each  of  these  matrices  is  of  the  form  (88)  and 

*2 

(89);  (compare  also  (10)  to  (lA)).  Accordingly,  the  submatrix  £  ,  given 

by 


*2 


'AA 

0 


0 

* 

'55 

0 


is  tiie  reciprocal  of  the  submatrix  S 


*2 


_  » 


0 

* 

^66  _ 
given  by 


(98) 


*2 


’AA 

0 


0 

* 

Ss 

0 


0 

0 

* 

*66 


(99) 


Since  these  matrices  are  diagonal,  their  reciprocity  is  ex¬ 
pressed  by  the  relations 


(a) 


29 


(b) 


* 


S 


55 


1 


> 


S*J  -  ^  .  (c)  ClOO) 

^66 

It  is  seen  that  when  (100)  is  introduced  into  (97),  lower 
bounds  for  C44,  C55,  and  are  obtained,  while  upper  bounds  are  given 
by  (96),  We  shall  now  adopt  the  physically  more  meaningful  notation 


*  * 

hi  ’  ^44  » 


*  * 

^’31  “  ^66  • 


(a) 

(b) 

(c)  (101) 


These  are  the  three  different  shear  moduli  o'  the  biaxially  reinforced  mat¬ 
erial.  The  subscripts  indicate  the  directions  of  shear.  Then  (96),  (97), 
and  (100)  provide  the  following  bounds  for  t’ e  shear  moduli: 


(a) 


(b) 


(c)  (102) 


It  will  be  recalled  that  the  extreme  right  and  left  sides  in 
the  inequalities  (102)  are  expressed  in  terms  of  the  single-layer  moduli 
and  compliances  and  the  reinforcement  angle  (|),  by  relations  of  type  (82). 
Therefore,  (102)  provides  lower  and  upper  bounds  for  the  effective  shear 
moduli  of  the  biaxially  stiffened  material.  Such  bounds  will  be  given 
In  detail  In  the  next  section. 

Unfortunately,  the  situation  for  the  other  six  effect^/e 
moduli  and  compliances  Is  not  so  simple,  since  the  matrices  (90)  and  (92) 
are  not  diagonal.  Furthermore,  the  effective  moduli  and  compliances 
entering  into  (90)  and  (92)  through  (87)  are  now\  contained  in  the  sub- 
matrices 


^-hl-  ^44 

’44 

^55 


r  ^  Si  <  S6 

’66 
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1 

'll 

^12 

S3 

* 

* 

* 

'12 

^22 

C23 

* 

* 

* 

13 

S3 

S3 

* 

* 

* 

11 

®12 

^13 

* 

* 

^22 

* 

12 

^23 

* 

* 

* 

13 

^23 

S33 

(103) 


(104) 


h  A 

Because  of  the  reciprocity  of  and  the  submatrices 

(103)  and  (104)  are  reciprocal.  However,  the  reciprocity  relations  are 
now  of  a  very  complicated  nature,  unlike  the  simple  relations  (100). 

The  mathematical  problem  with  which  we  are  now  confronted  is  as 
follows:  Given  that  the  matrices 

1  •!  *1 


•1  *1 

C  -  C  , 

-1  *1 

s  -  s  , 


(105) 

(106) 


are  positive  definite,  where  the  matrices  on  the  right  sides  are  of  the 
forms  (103)  and  (104),  and  furthermore  each  of  the  matrices  on  the  right 
sides  of  (105)  and  (106)  is  by  itself  positive  definite,  and  and  S 
are  reciprocal,  find  upper  and  lower  bounds  on  the  components  of  £*^. 

A  solution  to  such  a  problem  does  not  seem  to  have  appeared  in 
the  mathematical  literature.  The  present  bounding  method  will  be  used 
only  for  shear  moduli  bounding,  and  a  different  method  will  be  adopted  for 
the  other  moduli  and  compliances. 


8.  BOUiroS  FOR  SHEAR  MODULI  (LM  AND  CCM) 

The  right  sides  of  (102)  are  given  by  (82  h,i,j).  These  are 
now  rewritten  with  the  aid  of  (44)  to  (48).  From  (45),  (46),  and 
(82  h,i),  we  have 

C,  ,  ■  cos^ifi  G  +  sin^ij)  G^,  (107) 

44  a  t 

Ccc  “  sin^())  G  +  co8^()>  G  ,  (108) 

55  a  t 

From  (44)  to  (48)  and  (82  j),  we  have 

-  T  sln^2(t)  [E  +  G^  +  (1  -  2v  )^]  +  cos^2(})  G  .  (109) 

66  4  ^  a  t  t  a  a 

These  expressions  define  the  upper  bounds  in  (102). 

To  obtain  lower  bounds,  the  compliance  analogues  of  (82  h,i,j) 
are  needed.  But  here  great  caution  must  be  exercised  in  the  interpre¬ 
tation  of  the  six-by-six  matrix  notation.  Note  that  a  pure  shear  stress- 
strain  relation  for  orthotrcpic  material  is,  for  example, 


°12  "  ^^12  ^12  * 

which  becomes 

04  -  %  u  > 

in  the  notation  (81),  (85),  and  (86).  The  Inverse  is 

^4  "  ^44  *^4  » 

and  thus. 

However,  in  tensor  notation, 

^12  ■  ^^1212  °12  * 


(110) 

(111) 

(112) 

(113) 

(114) 


whereas  (112)  can  be  written  as 

^12  "  2  ^44  ^^12  • 


(115) 
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Thus,  from  (114)  and  (115), 


^44  “  ^^1212* 


(116) 


and  similarly, 


^55  “  ^^2323’ 


^66  “  ^^3131’ 


(117) 


in  contrast  to  relations  of  type  (81  b).  On  the  other  hand,  relations  of 
type  (81  a)  remain  of  the  same  form  for  compliances. 

In  velw  of  the  foregoing  comments  and  the  fact  that  (82)  Is 
really  fourth-rank  tensor  transformation,  (82  h,i,j)  assume  the  follow¬ 
ing  form  for  six-by-six  compliance  matrix  components: 


cos  4*  +  sin  ((i  ^^^22  ~  ^23^* 


’A4  "  ‘'44  '  "'“22 

.2.  ,  _ 2^  o/c*  _  Q'  ^ 

22  ^23^  ' 

-  4cos^<ti  sln^4i  +  S22  ■  281^2^ 


S55  =  s 


+  (sln^^i  -  cos^(J>)^  .  (j.) 


(118) 


The  compliances  in  (118)  are  defined  by  the  inversion  of  (42)  and  (43), 
which  may  be  written  as 


^ll  “ 

^ll  ^ll  +  ^12 

CM 
•  CM 

1  D 

^^l2 

°33  » 

(a) 

^22  - 

Su  "ll  ^ 

022 

^^23 

O33  , 

(b) 

633  - 

S|2  o[i  +  S23 

°22 

+  SJ2 

°33  * 

(c) 

2^2 

-s;4^i2  » 

(a) 

20^3 

■  — fs'  —  s'  ^ 
2^  22  23' 

0^3 

f 

(b) 

2e^3_ 

C  •  rt^ 

^44  °31  • 

(c) 
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Thus,  from  (120)  and  (43), 


.  ,  “  r.  ’ 


(121) 


c  '  —  S ' 

^22  23 


r'  r'  2G  • 
^22  "  ^23 


(122) 


By  Inoertlng  (121)  and  (122)  into  (118  a,b),  we  find 


-  cos^4>  ,  sin^4' 

S44  *  “g  ~g;~ 

a  t 


(123) 


2  2 

2  sin  41  .  cos  ^ 

^55  "  G  G  • 

a  t 


(124) 


By  combining  (123),  (124),  (107),  and  (108)  with  (102  a,b),  we  find  the 
following  pair  of  simple  bounds: 


r  2.  ,  2.>1 

cos  <P  .  sin  4 

G  G.  J  - 

La  t  J 


*  2  2 
G. ,  <  cos  41  G  +  sin  (}>  G  , 

•  ft  t 


(125) 


sin  4>  ,  COB 

G  G^ 

La  t 


♦  2  2 
2 —  I  <  G22  <  sin  4>  G^  +  cos  4)  G^  , 


(126) 


In  order  to  avoid  very  heavy  algebra,  a  different  procedure  will 
be  adopted  to  derive  an  expression  for  a  lower  bound  for  G^j.  For  this 
purpose,  (119)  will  be  rewritten  in  terms  of  Young's  moduli  and  Poisson's 
ratio  defined  in  Section  5.  Thus, 


o'  V  V 

”  1  1 1  ft  ^  I  ft  *“  f 

'^ll  “  E  "  E  '’22  ■  E  33  ' 
at  t 


-t  _a  — I  .  22  _  t 

*^22  “  "  E^  °11  Ej.  33  * 


+ 


(c) 


(127) 


It  is  seen  that  (127)  Is  merely  the  specialization  of  (13)  with 
(15)  and  (16)  to  a  transversely  Isotropic  material.  Comparing  (127)  with 
(119) ,  we  can  write 


sil 

1 

"  E  * 

a 

(a) 

S*2 

1 

(b) 

s* 

S2 

V 

a 

“  ”  ®t  * 

(c) 

s’ 

^23 

'-E 

■  ■  • 

(d) 

(128) 

Using  (121)  and  (128)  in  (118),  we  find 

S 


'66 


1  +  2v 

-  sin^2^  +  — I - ^ 

a  t 


^  1- 
+  cos  2^  —  . 

It 

a 


(129) 


By  combining  (109)  and  (129)  with  (102  c)  we  finally  obtain  bounds  on 
G31.  Thus 


|sln^2<ti 


1  +  2v 

_ I 


cos  2(^ 


-1 


<  G*J^  £  8in^2^  ^  [E^+G^  +  K^(l-2v^)^) 


+  cos  2^  G 


(130) 


Equations  (125),  (126),  and  (130)  thus  provide  lower  and  upper 
bounds  for  the  three  shear  moduli  of  a  biaxlally  fiber  reinforced  mat¬ 
erial  in  terms  of  the  transversely  isotropic  layer  moduli. 

Note  that  when  (ft  -  0,  the  blaxially  reinforced  material  degene¬ 
rates  into  a  unlaxlally  stiffened  material  with  fibers  in  the  x^  direction 
(Figure  2).  In  this  event,  (125),  (126),  and  (130)  reduce  to 


* 


Si  “  s  * 

as  they  should. 
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On  the  other  hand,  when  (j)  -  90®,  t’le  biaxially  reinforced 
material  degenerateo  into  a  uniaxially  stiffened  material  with  fibers  in 
the  X3  direction.  In  this  event,  (125),  (126),  and  (130)  reduce  to 


(132) 


as  they  should. 

We  now  consider  the  situation  when  the  layers  are  described  by 
the  composite  cylinder  model  (CCM) ,  Figure  3b.  The  bounding  method  to 
be  used  has  been  described  in  general  terms  in  Section  4.  Fortunately, 
it  is  not  necessary  to  give  any  further  analysis  based  on  this  method,  for 
it  turns  out  that  an  entirely  equivalent  procedure  is  to  take  for  elastic 
moduli  of  the  layers  those  predicted  by  Hashin  and  Rosen  [1],  and  to  use 
these  in  the  LM  bounding  method. 


The  truth  of  this  statement  follows  from  the  method  of  analysis 
employed  in  [1].  If  average  shear  strains  or  shear  stresses  are  applied 
to  a  uniaxial  fiber  composite  which  is  described  by  the  CCM,  then  the 
strain  or  stress  energies  stored  in  it  are  given  by  the  usual  expressions 
for  strain  or  stress  energies  in  homogeneous  materials,  with  the  effective 
shear  moduli  replacing  the  homogeneous  shear  modu'*!.  Therefore,  the  LM 
bounding  method  which  was  based  on  a  homogeneous  layer  becomes  immediately 
applicable  for  CCM  layers,  with  values  predicted  by  the  CCM  method  used  for 
the  layer  moduli. 


Turning  first  to  (125)  and  (126)  it  is  important  to  note  that 
the  CCM  method  gives  a  closed  expression  for  Gg  but  only  lower  and  upper 
bounds  for  G^.  Denoting  such  bounds  on  Gt  by  G^“^and  G^+) ,  we  note  the 
following  relations: 


Eco8^i(i  8in^(H~^  ^  I  cos^<|i 

A  G<-lJ  -L'-a 


t  - 


Eln^i^  co8^<(>  ^  <1  sin^iji 

G<-U  -L"a 


+  sln^4>  , 

(a) 

8in^(^”|■^ 

^’t  J 

(b) 

+  C08^(^  , 

(c) 

COS^^"^ 

^t  J  • 

(d) 

(133) 
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The  best  bounds  available  for  of  the  CCM  are  given  by  (64) 
and  (65).  Thus,  with  the  use  of  these  in  (133),  the  bounds  (125)  and 
(126)  assume  the  form 


fcos^^  ^  p* 

b. 


2  2  (+) 
cos  ij)  G^  +  sin  ^  ®t(c)  * 


2  2  (+) 
sin  4)  G  +  cos  4’  G  >  »  • 

a  t('') 


(134) 


(135) 


Similar  remarks  concern  (130),  G^  appears  again  on  the  right 
side,  while  all  other  moduli  on  this  side  have  closed  expressions  on  the 
basis  of  the  CCM.  On  the  left  side,  all  the  moduli  except  Et  have  closed 
expressions.  Bounds  on  can  be  expressed  in  terms  of  bounds  on  G^  by 
(49),  Note  in  this  respect  (62),  Let 


(136) 


where  the  functional  relationship  in  (136)  is  given  by  (49)  with  in¬ 
equality  arguments  similar  to  those  used  in  the  derivation  of  (134)  and 
(135),  We  find  on  the  basis  of  (130), 


8ln^2(^ 


/  'a  J  - 


<  i  8ln22*pE^  +  +  Kjd  -  2v^)  J  +  coa22*  .  (137) 


I 


9.  NUMERICAL  RESULTS  FOr>  SHEAR  MODULI  BOUNDS 
By  utilizing  the  theoretical  results  derived  In  Section  8, 


bounds  on  the  effective  shear  moduli 
Expressions  for  bounds  are 


will  be  calculated  In  the  following, 
given  by  (13A),  (135),  and  (137); 

G*  <  cos^(J)  G  +  sln^({)  Gft\  (134) 

12  -  ^  a  ^  t(c) 

G23  <  sin  <}>  +  cos^(|)  (135) 


^  2^  .  ,1  .  l+2va  ^  .  cos^24' 

sin  2MeJ  +  -^)) 

“  E,.  a 

-  t  _ 


r* 

i  Si 


1  2  2  2 
<  i  sin  2<t  [E^  +  gJ^^j  +  Kj.  (l-2Va)  ]  +  cos  2^  G^ 


(137) 


The  moduli  appearing  In  these  equations  are  calculated  from 


expressions  given  In  Section  5. 


E  -  E  V  +  E^v 
a  mm  f 


(62) 
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-G  r  1  ~  V  aO  . 

’t(c)  %L  1  -  2v  \  J  ’ 

^  m 


(138) 


In  (137),  is  obtained  from  (49  a),  replacing  by  G^”^, 


Kj  +  ♦  G*"’ 


(139) 


'I' 


1  + 


4K 

t  a 


(49  b) 


The  elastic  constants  both  for  the  fiber  and  the  matrix  materials  (de¬ 
noted  by  the  subscripts  f  and  m,  respectively)  are  given  by 


E 

2(1  +  v)  ’ 

(a) 

vE 

(b) 

(1  +  v) (l-2v)  » 

A  +  G  . 

(c) 

(140) 

.(+), 


.-c 


An  expression  for  Gt(c)  is  given  in  [1],  A,'"  can  be  obtained 

from  the  solution  of  the  following  six  linear  algebraic  equations,  [1]: 
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m 

-n 
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^5 
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1 
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0 

n 

TE 

0 

3-2v 

m 
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m 

3-2v^ 

^6 

(lAl) 


where 


G,/G 
f '  m 


(142) 


The  results  given  in  this  section  have  been  programmed  for 
the  IBM  709A  digital  computer.  For  composite  materials  whose  properties 
are  listed  below,  sample  calculations  have  been  performed. 


Fiber  glass/Epoxy 

-  10.5  X  10^  psi 

-  0.2 


E  ■  0.5  X  10^  psi 
m 

V  -  0.3  X  10^ 
m 


Boron/Magnesium 
E^  ■  60  X  10^  psi 


E  ■  6.5  X  10^  psi 
m 


V 


f 


0.3 


V  -  0.3 

n 
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The  results  are  shown  In  graphical  form  In  Figures  4  through 
39.  Figures  4  through  21  are  for  fiber  glass/epoxy;  and  Figures  22 
through  39  for  boron/magneslum.  Shown  are  variations  of  shear  moduli 
with  angle  of  reinforcement  (fi  for  fixed  volume  fraction  and  variations 
with  fiber  volume  fractions  for  fixed  angle  of  reinforcement  4*. 

h  it 

It  is  seen  that  the  bounds  on  G^2  ^23  extremely  close, 

while  those  on  G*^  are  further  apart. 
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SHEAR  MODULUS  (10**  psi) 


* 

Figure  4.  Shear  Modulus  VS.  Volume  Fraction 
(|)  *  30  degrees 
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SHEAR  MODULUS  (10®  psi ) 


VOLUME  FRACTION 


Figure  5.  Shear  Modulus  G23  vs.  Volume  Fraction 
41  =  30  degrees 
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SHEAR  MODULUS  (10®  psi) 


* 

Figure  6.  Shear  Modulus  ^12  vs.  Volume  Fraction 
4,  =  30  degrees 
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VOLUME  FRACTION 


•k 

Figure  7,  Shear  Modulus  Gg-j  vs.  Volume  Fraction 
4,  =  45  degrees 
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SHEAR  HOOULUS  (lO**  psi) 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

VOLUME  FRACTION 

★ 

Figure  8.  Shear  Modulus  ^23  vs.  Volume  Fraction 
(fi  =  45  degrees 
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SHEAR  MODULUS  (10®  psi) 


* 

Figure  9,  Shear  Modulus  G12  vs.  Volume  Fraction 
({i  =  45  degrees 
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SHEAR  HODULUS  (10°  psi) 


f 

f  • 


* 

Figure  10.  Shear  Modulus  6^^  vs.  Reinforcement  Angle 
Volume  Fraction  =  0.3 
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SHEAR  MODULUS  (10°  psi) 


Figure  12, 


* 

Shear  Modulus  G^2  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,3 
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REINFORCEMENT  AN6LE.  ♦  (deg) 


K 

Figure  13.  Shear  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =0.5 


REINFORCEMENT  ANGLE.  ♦  (deg) 


« 

Figure  15,  Shear  Modulus  6^2  ^s.  Reinforcement  Angle 
Volume  Fraction  =  0.5 


SHEAR  MODULUS  (10®  psi) 


fiber  glass/epoxy 

X  >  upper  bound 
0  »  lower  bound 


REINFORCEMENT  ANGLE.  4  (deg) 


Figure  17,  Shear  Modulus  G23  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,7 


SHEAR  MODULUS  (10®  psi  ) 
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0.4 
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* 

Figure  18,  Shear  Modulus  G^2  vs.  Reinforcement  Angle 
Volume  Fraction  =  0.7 


56 


REINFORCEMENT  ANGLE.  4  (deg) 


★ 

Figure  19.  Shear  Modulus  G^i  vs.  Reinforcement  Angle 
Volume  Fraction  =0.9 
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SHEAR  MOOUi-US  (10®  psi) 


REINFORCEMENT  ANGLE,  ♦  (deg) 


Figure  20,  Shear  Modulus  G23  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,9 
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20  40  60 

REINFORCEMENT  ANGLE.  4  (deg) 

* 

Figure  21,  Shear  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,9 


VOLUME  FRACTION 


Figure  22, 


* 

Shear  Modulus  G^-j  vs.  Volume  Fraction 
(|)  =  30  degrees 


60 


SHEAK  modulus  (10°  psi) 


I 


Figure  23. 


* 

Shear  Modulus  G23  vs.  Volume  Fraction 
(()  =  30  degrees 


61 


t 


I 

I 


1 


VOLUME  FRACTION 


* 

Figure  24.  Shear  Modulus  6^2  vs.  Volume  Fraction 
(j)  =  30  degrees 
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★ 

Figure  25,  Shear  Modulus  G^i  vs.  Volume  Fraction 
({>  =  45  degrees 
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SHEAR  MODULUS  (10®  psi) 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 
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Figure  26. 


* 

Shear  Modulus  G23  vs.  Volume  Fraction 
(|)  =  45  degrees 
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SHEAR  MODULUS  (10“  psi) 


20  40  60 

REINFORCEMENT  ANGLE,  ♦  (deg) 


* 

Figure  28,  Shear  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0.3 


* 

Figure  29,  Shear  Modulus  G23  vs.  Reinforcement  Angle 
Volume  Fraction  =0,3 
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20 


60 


Figure  30, 


40 

REINFORCEMENT  ANGLE,  «  (deg) 

★ 

Shear  Modulus  vs.  Reinforcement  Angl 
Volume  Fraction  =  0,3 


10°  psi) 
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40  60 

REINFORCEMENT  ANGLE,  ♦  (deg) 
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S’ 


Figure  33 


* 

.  Shear  Modulus 
Volume  Fraction  « 


vs.  Reinforcement  Angle 
0.5 
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Figure  34, 


* 

Shear  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,7 
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REINFORCEMENT  ANGLE. ^  (deg) 


Figure  35 


* 

,  Shear  Modulus  G23  vs.  Reinforcement  Angle 
Volume  Fraction  =  0.7 
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Figure  36. 


* 

Shear  Modulus  ^12  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,7 
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SHEAR  MODULUS  (10°  psi) 


0  20  40  60  80 

REINFORCEMENT  ANGLE,  ♦  (deg) 


Figure  38, 


★ 

Shear  Modulus  G23  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,9 
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Figure  39, 


★ 

Shear  Modulus  ''5*  Reinforcement  Angle 
Volume  Fraction  =  0,9 
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10.  COMPUTATION  OF  NORMAL  EFFECTIVE  ELASTIC  MODULI  FOR  LAYER  MODEL 


The  normal  effective  elastic  moduli  of  the  composite  are  here 
defined  as  the  six  moduli  C*j  =  Cj^,  1,  j  ■  1,2,3  which  relate  the  normal 
average  stresses  to  the  normal  average  strains;  see  (10),  It  will  be 
shown  here  that  if  the  composite  is  described  by  the  layer  model  with 
transversely  isotropic  layers,  then  the  normal  Cjj  can  be  exactly  com¬ 
puted  in  terms  of  the  layer  moduli. 


Let  a  specimen  of  the  composite,  as 
subjected  to  the  boundary  disp? acements 

shown  in 

Figure  1,  be 

Ui(S)  « 

^2(3)  “  ^22  ^2* 

(^) 

(b) 

U3(S)  -  C33  X3 , 

(c) 

(143) 

where  eii,  G22, 
result  (3)  tMt 

and  ^23  constant*  It  then 

the  average  strains  are 

follows 

from  the  general 

hi  “  ^11’ 

(a) 

h2  “  ^^22* 

(b) 

h3  “  h3  * 

(c) 

h2  "  h3  "  *^31  “  ° 

(d) 

(144) 

The  average  stresses  are  then  related  to  the  average  strains  by  (10) .  We 
shall  show  here  that  in  the  present  case  the  average  stresses  can  actually 
be  computed.  These  then  lead  via  (10)  to  expressions  for  the  C*.  appear¬ 
ing  in  (10), 

We  look  for  an  exact  elasticity  solution  of  the  layered  material 
under  boundary  conditions  (143;,  We  make  the  guess  that  the  internal 
displacements  throughout  all  layers  are  of  the  form  (143).  Thus, 


Uj^(x)  -  Xj^, 

(a) 

^2(21)  “  ^22  ^2* 

(b) 

u^()c)  “  ^33  ^3* 

(c) 

(145) 
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For  (145)  to  be  an  exact  elasticity  solution  of  the  layered 
material,  it  must  satisfy  the  following  conditions: 

(1)  The  boundary  conditions  (143), 

(2)  The  differential  equations  of  elasticity  theory 
throughout  all  layers, 

(3)  Displacement  and  traction  continuity  at  layer  interfaces. 

Now  (143)  is  obviously  satisfied,  since  the  choice  (145)  was 
made  for  that  purpose. 

The  differential  equations  of  elasticity  of  the  layers  are 
second-order  partial  differential  equations.  Since  the  equations  (145) 
are  linear,  they  satisfy  these  equations  trivially. 

Displacement  continuity  at  layer  interfaces  is  obviously  sat¬ 
isfied,  since  (145)  is  by  choice  the  same  continuous  field  in  all  the 
layers.  Thus,  only  traction  continuity  at  layer  interfaces  remains  to 
be  investigated.  Since  the  normal  to  the  layer  interface  is  in  the  X2 
direction.  Figure  1,  the  components  of  the  traction  at  the  interface  are 
°2l»  ®22»  ®23*  Thus,  we  must  have 


21 


'22 


23 


)  continuous  at  layer  interfaces 


(146) 


This  will  now  be  investigated. 


First,  compute  the  strains  and  In  the  layers.  The 


strains  associated  with  (145)  are 


ij 


’ll 

0  e 

0 


22 

0  E 


0 

0 

33 


(147) 


Using  the  transfoimatlons  (70)  and  (72),  we  find 
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ai 

2 

cos  (fi  + 

G33 

sin^4* 

■'a 

9 

(a) 

^22 

- 

^22 

-  eJJ 

9 

(b) 

£  * 

3; 

m 

2 

sin  <(i  + 

633 

cos  <(| 

-  e" 

33 

9 

(c) 

£  * 

1:. 

- 

(‘33  •  'll) 

sin 

41  cos  4 

-  -e" 

13 

9 

(d) 

0 

.  e*' 

12 

9 

(e) 

*^23 

- 

0 

-  e" 

23 

• 

(f) 

(148) 


The  stresses  in  the  primed  layers,  Figure  2a,  are  given  by 
(42).  Similarly,  the  stresses  in  the  double  primed  layers.  Figure  2b, 
are  given  by  (42)  with  replaced  by  eV 


11 


22 


33 


12 


23 


13 


Ij* 

I  • 
11* 

I I 

22* 

O33, 


j” 

12 


t  •  • 


■  -0 


23 

1 1 


Using  (148),  we  find 
(a) 


0. 

0, 


13* 


(b) 

(c) 

(d) 

(e) 

(f) 


(150) 


In  order  to  investigate  the  continuity  conditions  (146),  it  is 
necessary  to  transform  the  layer  stresses  (149)  back  to  the  composite 
coordinate  system  Xj^,  X2,  X3.  For  this  purpose,  we  use  the  transforma¬ 
tions  (71)  and  (73),  As  has  been  done  previously,  we  denote  the  left 
sides  of  (71)  by  and  the  left  sides  of  (73)  by  •  This  notation 

means  stresses  in  layers  referred  to  the  x^,  x^,  x^  system  of  axes  ex¬ 
pressed  in  terms  of  stress  components  referred  to  the  xj^,  x^,  X3  or  x'^, 
x^  systems,  respectively. 


We  now  have,  in  view  of  (149)^ 

''^21-  ”'^21  - 

(a) 

'‘^22  "  *^22  “  °22  “  "°22* 

(b) 

’'^23  “  ”^23  - 

(c) 

(150) 
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Therefore,  the  continuity  conditions  (146)  are  satisfied.  We  thus  con¬ 
clude  that  the  displacement  field  (145)  is  an  exact  elasticity  solution 
for  the  layered  composite. 


% 


We  shall  now  need  the  average  stresses  <^22’  ®33* 

this  purpose,  we  first  compute  the  uniform  layer  stresses  ’o22» 

'032  and  "^22*  "°33*  a,b,c)  and  (149),  we  find 


11 


22 


2  2 

=  cos  +  0^^  sin  (})  -  2oj^  sin  1(1  cos  ii,  (g) 


22’ 


(b) 


2  2 

=  o'  sin  i{i  +  o'  cos  ())  +  2o'  sin  <\i  cos  ,  (c) 


'33  ”11  "  ■  ”33  13 

Similarly,  from  (73  a,b,d)  and  (149),  we  have 


(151) 


'11 


22 


33 


'11 

j' ' 

22’ 

t  •  • 


cos  (fi  +  0^2  sin  4>  +  sin  cos  41,  (a) 


13 


(b) 


2  2 

o||  sin  4  +  0^2  cos  4>  -  2o|j  sin  4  ^os  4>,  (c)  (152) 


In  view  of  (149),  it  is  seen  from  (151)  and  (152)  that 


^  n  ■  ^  *  0 

°n  11’ 


'22 


1 1, 


'22' 


33 


33' 


(a) 

(b) 

(c)  (153) 


where  (153b)  has  already  been  found  in  connection  with  the  continuity 
conditions. 


Since  the  stresses  (153)  are  uniform  and  since  both  kinds  of 
layers  are  of  equal  volumes,  we  have 


°11  “  '°11’ 

(a) 

022  ■  '022’ 

(b) 

°33  “  '°33* 

(O 

(154) 

We  first  compute  the  layer  stresses  in  terms  of  tii , 

do  this,  (148)  is  used  in  (42),  yielding 

£22*  £33. 

To 

I 

m 


1 


cos’ip  +  CJ2  sin^(j)) 


O  0 

+  (c|^  sin  41  +  0^2  cos  4)  ,  (a) 


0*2  =  (c;2  cos  4  +  ^3  sin  4>)  +  0*2  ^22 

+  (CJ2  sin^4  +  C^3  cos^(})]  C33  ,  (b) 

oi_  *  ^^12  "*"  ^22  ^11  ^23  ^22 


'33  '12 


2  2  ^ 

+  (C^2  ***  ■'■  *^22  ^33  » 


°l3  “  ^^33  "  ^11^  ^ 


(d)  (155) 


Now  (155)  Is  subat.ituted  into  (151),  which  in  view  of  (15A)  are  also 
the  average  stresses.  After  rearrangement,  we  have 


4  A  2  2 

cos  ^  +  C^2  si*'  'P  +  ^^12 


+  8ln^(|)  cos^i^J  £33  +  ^C|2  cos^(()  sln^(t)Je22 

+  +  C22)  sln^*  008^4  +  C|2  (sin^^i  +  cos^4i) 

-  ACJ^  sln^4  cos^^J  , 


2  2 

-  AC,’,  sin  *  cos  4 
AA 


°22  •  ^2 


+  C^3  sin^^J 


^11  ■'■  ^22  ^22 


+  jc32  sin^4  +  co8^i|> 


^33* 


3  "  ^^^il  '*'  ^22^  sin^cfi  co8^((i  +  0^2  (sin^((i  +  cos^^) 

-  AC^^  sin^^i  cos^i^  +  ^C|2  sin^(ji  +  €^3  co8^(|i^ 

■i  jc|3  sln^4i  +  C^2  ^08^4  +  2CJ2  cos^4>  sin^ 

+  AC^^  cos^ij)  sin^^i^  e33»  (c)  (156) 


82 


Comparing  (156)  with  (10),  we  see  that  the  effective  elastic  moduli  are 


/  A  2  2 

C*^  -  cos^4i  +  C'2  sin  %  +  '1“  ♦ 

2  2 

+  cos  (J)  sin  (J), 


C*2  •=  C*^  -  CJ2  cos%  +  C’3  sin^^. 


(157) 

(158) 


-  4C!.  co8^(j)  sin  ((>, 


C22  "  Cj2* 

C^3  -  C*2  -  C|2  8in^(^  +  cos^*, 

i  i  2  2 

^22  ^  ^^12 

2  2 

+  4CJ,  cos  (J)  sin  41, 


(159) 

(160) 
(161) 


(162) 


If  the  elastic  moduli  of  the  layers  are  knovm  from  experiment 
or  theory,  then  (157)  through  (162)  enable  one  to  compute  the  six  effec¬ 
tive  elastic  moduli  in  (10).  However,  experimental  determination  of  the 
C].  is  difficult.  Therefore,  it  is  best  to  proceed  on  the  basis  of 
theoretical  knowledge  about  the  It  will  be  shown  in  the  next  section 

how  (157)  through  (162)  can  be  used  to  bound  the  moduli  C*  . 


11.  BOUNDS  FOR  NORMAL  EFFECTIVE  ELASTIC 
MODULI  FOR  COMPOSITE  CYLINDER  MODEL 


At  the  present  time,  no  geometrical  model  of  a  uniaxially  fiber 
reinforced  material  for  which  all  five  elastic  moduli  can  be  rigorously 
determined  is  available. 


The  composite  cylinder  assemblage  model  of  Hashln  and  Rosen  [1] 
comes  closest  to  such  a  model,  since  it  permits  rigorous  computation  in 
closed  form  of  four  of  the  five  elastic  moduli;  see  Section  5.  Therefore, 
we  shall  assume  that  the  layers  are  described  by  this  model,  and  we  shall 
now  proceed  to  exploit  the  theory  developed  in  Section  10  for  this  case. 


First,  the  layer  moduli  appearing  in  (157)  through  (162)  are 
rewritten  by  use  of  (AA)  through  (A8).  Thus, 


c' 

E  + 

A  K 

11 

a 

t 

c’ 

2  K 

V 

12 

t 

a  ♦ 

c* 

K  + 

G  . 

22 

t 

t  » 

C' 

m 

K  - 

23 

t 

t  » 

c* 

^AA 

- 

"a  • 

(163) 

(16A) 

(165) 

(166) 

(167) 


Inserting  these  expressions  into  (157)  through  (162)^  we  have 

cjj^  ■  E^  cos^(^  +  8in^(;))^  (168) 

2  2  A 

+  A  cos  <p  sin  (Ji  +  sin  (P  ,(169) 
C*2  ■  (2Vg  cos^<|)  +  sin^4i)  -  G^  sin^((i  ,  (170) 

C*2  ■  Eg  cos^^i  sin^<;i  +  Kj(2v^  cos^<|>  +  sin^<(i)  (2Vg  8in^(}i  +  cos^(fi) 

-  A  Gg  co8^(P  sin^<p  +  G^.  co8^((i  sin^(>  ,  (171) 

0^2  -  Kt  +  Gj.  ,  (172) 

^23  “  ^  co8^(^)  -  G^  co8^(<i  ,  (173) 


8A 


(17A) 


C33  “  Eg  sin^(J>  +  cos^(j>)^  ^  ^a  8in^()) 

+  G^  cos^(fi  . 

The  following  problem  now  arises:  While  the  theory  which  has  been  given 
in  [1]  gives  rigorous  results  in  closed  form  for  the  moduli  Eg,  Vg,  K  , 
and  Ga»  only  lower  and  upper  bounds  for  Gf.  are  available;  see  (6AJ  and 
(65),  It  is  seen  that  G^.  appears  in  every  one  of  the  Equations  (168) 
through  (17A).  Therefore,  only  the  upper  and  lower  bounds  can  be  calcu¬ 
lated  for  the  moduli  in  these  equations. 

A  rigorous  procedure  for  bound  construction  for  the  effective 
moduli  appearing  on  the  left  sides  of  (168)  through  (174)  would  consist 
of  the  following  steps.  Displacements  of  the  form  of  (143)  are  applied 
as  boundary  conditions  to  the  surfaces  of  all  composite  cylinders  of  which 
the  layers  are  assumed  to  consist,  see  Figure  3.  A  solution  can  be  ob¬ 
tained  for  a  typical  composite  cylinder  and,  therefore,  for  each  of  the 
cylinders  in  the  layers.  The  assemblage  of  displacement  fields  in  all 
the  cylinders  then  becomes  an  admlssable  field  for  the  principal  of  mini¬ 
mum  potential  energy. 

Then  the  sum  of  the  strain  energies  in  all  cylinders  subjected 

to  (143)  on  their  surfaces  is  larger  than  the  composite  strain  energy, 

and  an  inequality  of  the  general  type  (35)  is  obtained.  The  right-hand 
side  of  (35)  in  the  present  case  is  a  quadratic  form  in  which  the  vari¬ 
ables  are  £22*  £33  appearing  in  (143),  and  the  coefficients  are  ex¬ 
pressed  in  terms  of  the  CCM  moduli  Eg,  Kj.,  Vg,  Gg,  and  which  is  given 

by  (65),  The  left-hand  side  of  (35)  now  contains  the  unknown  normal  ef¬ 
fective  moduli  of  the  composite  written  on  the  left-hand  sides  of  (157) 
through  (162). 

An  entirely  similar  procedure  can  be  used  with  traction  boundary 
conditions  involving  only  normal  stresses  ,  022*  and  033  with  respect 
to  the  composite  system  of  axes.  The  general  procedure  has  been  described 
in  Section  4  and  now  leads  to  an  inequality  of  type  (41),  The  inequalities 

of  type  (35)^and  (41)  have  to  be  exploited  to  find  the  upper  and  lower 

bounds  for  .  As  has  been  mentioned  in  Sections  4  and  7,  it  is  at 
present  not  known  how  to  solve  such  inequalities.  Therefore,  we  shall 
have  to  adopt  a  different  procedure. 

It  will  simply  be  assumed  that  if  the  layers  are  described  by 
the  CCM,  we  can  replace  the  whole  composite  by  one  whose  layers  are  homo¬ 
geneous  and  have  the  uniaxial  CCM  elastic  moduli.  It  is  quite  clear  that 
the  changes  thus  Introduced  into  the  effective  elastic  moduli  of  the  CCM 
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biaxial  composite  will  be  Insignificant,  If  each  layer  contains  an  ap¬ 
preciable  number  of  fibers  through  Its  thickness.  On  the  basis  of  this 
assumption,  bounds  for  C*j  are  obtained  immediately  by  replacement  of 
in  (168)  through  (174)  by  its  lower  or  upper  bounds,  as  given  by  (64) 
and  (65).  It  is  seen  that  in  Equations  (168),  (170),  (171) ^ and  (174), 
appears  with  positive  coefficients.  Therefore,  these  expressions  are 
a  monotonlcally  increasing  function  of  G^.  Consequently,  insertion  of 
upper  bound  values  for  G^  results  in  upper  bounds  for  the  left  sides, 
and  insertion  of  lower  bound  values  results  in  lower  bounds  for  the 
left  sides.  In  equations  (169)  and  (173)  Gj.  appears  with  negative  coef¬ 
ficients,  so  these  expressions  monotonlcally  decrease  with  increasing  G^., 
Therefore,  upper  bound  values  of  G^  give  lower  bounds  for  the  left  sides, 
and  lower  bound  values  of  Gj.  give  upper  bounds  for  the  left  sides. 


As  a  result  of  these  considerations,  we  can  now  write  the 
following  set  of  bounds  for  the  effective  moduli  0^4  on  the  left  side  of 
(168)  through  (174). 


'(+) 


-11 


<  p*  p*(+) 

^11  f  ^11  i  ^ll 

(a) 

(2va  cos^(j)  +  sin^((()^  +  4  G 

2  2 
a  cos  (ji  sin 

+ 

G^~^  sin^(ji  , 

(b) 

(2va  cos^((>  +  8ln^(ji)^ 

^  cos^((i  sin^cfi  +  sin^^i  , 

(c) 

cl*-)  i  i  ci^) , 

(a) 

(2va  co8^(j»  +  sin^ifi)  - 

sin^ij)  , 

(b) 

(2va  cos^iji  +  sln^ij))  -  G^~^ 

2 

sin  1^  , 

(c) 

^13  -  ^13  -  ^13  ’ 

(a) 

,*(+)  . 
-13 


cos^(J)  sin^<ti  +  (2Vg  cos^())  +  sin^(<i)(2va  sin^())  +  cos^(|>) 

2 

^  roj 

a 


-  4  G_  cos^(|i  8ln^(}»  +  cob^<P  sin^^i  ,  (b) 


E^  cos^({>  sin^(j)  +  (2Va  cos^<p  +  8in^(|i)(2v^  8in^(ji  +  cos^(j)) 

-  4  G  cos^fli  8in^((i  +  cos^4)  8in^(j)  ,  (c)  (177) 

a  ^ 
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^22  -  ^22  -  ^22  ’ 

(a) 

^22  “  ^t  ^  * 

(b) 

^22  ~  h  ’ 

(c) 

(178) 

r*(-)  ,  *  .  p*(+> 

23  -  23  -  23  ’ 

(a) 

€*3"^  =  Kt  (2v^  sln^4>  +  cos^ij))  - 

2 

cos  4>  , 

(b) 

€23^^  “  Kj  cos^ifi)  “  ^ 

cos  ij)  , 

(c) 

(179) 

S3  -  33  -  33 

f 

(a) 

sin^(<>  +  Kj.  (2V3  Sin^4>  +  cos^tp)^  ^  cos^p  Bin^(p 

+  C08^(|>  ,  (b) 

-  E  8in^((i+  K,  (2v-  sin^  <j>  +  cos^4')^  +  4  G  cos^'t'  sin^ip 
Ji  a  t  “  a 

+  G^"*"^  cos^iji  .  (c)  (180) 

The  bounds  (175)  through  (180)  have  been  numerically  evaluated 
for  fiber  glass/epoxy  and  boron/magnesium  composites  whose  constituent 
elastic  properties  have  been  listed  in  Section  9.  The  uniaxial  layer 
composite  moduli  have  again  been  computed  on  the  basis  of  the  results 
listed  in  Section  9, 

Bounds  for  fiber  glass/epoxy  composites  for  varying  angles  of 
reinforcement  and  fiber  volume  fractions  are  given  in  Figures  40  through 
57|  while  similar  bounds  for  boron/magnesium  composites  are  given  in 
Figures  58  through  72. 

It  is  seen  that  in  a  considerable  number  of  cases,  the  bounds 
are  so  close  that  they  are  practically  indistinguishable. 
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Figure  40, 


* 

Effective  Modulus  C^-j  vs.  Volume  Fraction 
(|)  =  45  degrees 
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VOLUME  FRACTION 


Figure  41 . 


* 

Effective  Modulus  C22  vs.  Volume  Fraction 
♦  =  45  degrees 
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SHEAR  MODULUS  (10®  psi) 
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SHEAR  MODULUS  (10°  psi) 


* 

Figure  43,  Effective  Modulus  Cii  vs.  Reinforcement  Angle 
Volume  Fraction  *  0,2 
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SHEAR  MODULUS  {10“  psi ) 
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Figure  44, 


* 

Effective  Modulus  0^3  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,2 
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SHEAR  MODULUS  (TO®  psi) 
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Figure  45.  Effective  Modulus  C22  vs.  Reinforcement  Angle 
Volume  Fraction  =0,2 
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SHEAR  MODULUS  (TO®  psi) 
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REINFORCEMENT  ANGLE,  ♦  (deg) 

* 

Figure  46,  Effective  Modulus  C-ji  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,5 
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* 

Figure  47,  Effective  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0.5 
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SHEAR  MODULUS  (10®  psi) 


* 

Figure  48,  Effective  Modulus  C22  vs.  Reinforcement  Angle 
Volume  Fraction  =0,5 
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SHEAR  MODULUS  (10®  psi) 


VOLUME  FRACTION 


* 

Figure  49.  Effective  Modulus  C-jg  vs.  Volume  Fraction 
(j,  =  45  degrees 
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SHEAR  MODULUS  (10°  psi) 
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Figure  50, 

* 

Effective  Modulus  Cp  vs.  Volume 

Fraction 

4)  =  45  degrees 
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SHEAR  MODULUS  (10*'  psi) 


* 

Figure  51,  Effective  Modulus  C23  vs.  Volume  Fraction 
(^  =  45  degrees 
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SHEAR  MODULUS  (10“  psi) 


REINFORCEMENT  ANGLE,  ^  (deg) 

★ 

Figure  52.  Effective  Modulus  C-jg  vs.  Reinforcement  Angle 
Volume  Fraction  =0.2 
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* 

Figure  53,  Effective  Modulus  vs.  Reinforcement  Angle 
Volime  Fraction  =  0.2 
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★ 

Figure  54,  Effective  Modulus  C23  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,2 
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Figure  55 


REINFORCEMENT  ANGLE,  ♦  (deg) 


.  Effective  Modulus  C*2  vs.  Reinforcement  Angle 
Volume  Fraction  0,5 
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Figure  56, 


Effective  Modulus  C23  vs.  Reinforcement  Angle 
Volume  Fraction  =  0.5 
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SHEAR  MODULUS  (10°  psi) 


REINFORCEMENT  ANGLE,  ♦  (deg) 


★ 

Figure  57,  Effective  Modulus  C-jj  vs.  Reinforcement  Angle 
Volume  Fraction  =0,5 
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SHEAR  MODULUS  (10®  psi) 
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SHEAR  MODULUS  (10^  psi) 
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Figure  59. 

* 

Effective  Modulus  C22  vs.  Volume  Fraction 

41  =  45  degrees 
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* 

Figure  60,  Effective  Modulus  C33  vs.  Volume  Fraction 
({,  =  45  degrees 
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Figure  62,  Effective  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,2 
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Figure  63.  Effective  Modulus  C2Z  Reinforcement  Angle 
Volume  Fraction  *  0,2 
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Figure  64.  Effective  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,5 
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Figure  65,  Effective  Modulus  vs.  Reinforcement  Angle 
Volume  Fraction  =  0,5 
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Figure  66.  Effective  Modulus  ''5*  Reinforcement  Angle 
Volume  Fraction  =  0.5 
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Figure  68,  Effective  Modulus  ^13  vs.  Volume  Fraction 
<|)  =  45  degrees 
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VOLUME  FRACTION 
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Figure  69,  Effective  Modulus  C23  vs.  Volume  Fraction 
4,  =  45  degrees 
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SHEAR  MODULUS  (10  pSi) 
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Figure  71,  Effective  Modulus  C23  vs.  Reinforcement  Angle 
Volume  Fraction  *  0,2 
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* 

Figure  72.  Effective  Modulus  C-jj  vs.  Reinforcement  Angle 
Volume  Fraction  *0,2 
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12.  BOUNDS  FOR  NORMAL  EFFECTIVE  EIASTIC 

COMPLIANCES  FOR  COMPOSITE  CYLINDER  MODEL 

The  problem  now  arises:  Can  bounds  on  the  normal  effective 
compliances  be  generated  on  the  basis  of  the  knowti  bounds  on  the  effec¬ 
tive  elastic  moduli? 


Analytical  treatment  of  this  problem  Is  unfortunately  prohib¬ 
itively  difficult.  The  analytical  procedure  would  be  to  express  the  com- 

iHt 

pllance  S^j  by  Inversion  of  the  matrix  C^j  and  then  to  Investigate  whether 

the  resulting  expressions  are  monotonlcally  Increasing  or  decreasing  In 

A  * 

terms  of  C , .  variations.  The  complexity  of  the  expressions  for  S,.  In 

terms  of  C^^  rules  out  such  an  Investigation. 

The  procedure  which  Is  here  adopted  Is  a  numerical  one.  The 

compliances  C^j  as  given  by  (168)  through  (174)  are  computed  for  a  boron/ 

epoxy  composite.  The  elastic  properties  of  matrix  and  fibers  are  given  by 

E  -  0.5  x  10^  psla  E,  »  57.5  x  10^  psl 

m  I 

V  -  0.35  V,  -  0.2 

m  t 

The  calculations  will  be  performed  for  four  values  of  G^.  These 

are  the  lower  and  upper  bounds  on  G^(l.e.,  ^  and  G^"**^),  as  given  by 

unlaxlally  composite  theory,  as  described  above,  and  by  two  Intermediate 

values 


.  g/->  +  i 

-  <=/->) . 

(181) 

-  g/->)  . 

(182) 

As  has  been  pointed  out  In  Section  11,  (168),  (170),  (171),  and 
(174)  are  monotonlcally  Increasing  functions  of  G^,  while  (169)  and  U73) 

are  monotonlcally  decreasing  functions  of  G^.  Thus,  for  the  values  G^“^ 

(+)  *  ^  ^ 
and  G'  ,  the  bounds  for  the  C. .  are  obtained,  while  the  values  (181)  and 

(182)  give  Intermediate  values.  C  results  for  G^  '  are  shown  in  Table  l; 

(1)  (2)  ^ 
for  G^  In  Table  IIIj  for  gJ  in  Table  IV;  and  for  G^  ,  in  Table  II. 

A 

Now  the  matrix  Is  numerically  Inverted  In  each  case,  whereby 
four  sets  of  values  are  obtained.  These  results  are  shown  in  Tables 
V  through  VIII. 
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TV 

According  to  the  numerical  results,  the  are  monotonlc 
functions  of  G^.  It  Is  seen  that  lower /upper  bounds  for  C, ,  generate 
upper/lower  bounds  for  the  diagonal  compliances  ^22*  ^33* 

cause  of  (15)  upper /lower  bounds  on  S-«,  and  S--  correspond  to 
lower/upper  bounds  on  Ej^,  E2,  and  E^.  Thus,  lower/upper  bounds  for 
generate  lower/upper  bounds  for  E^^,  E2,  and  E^. 

*  *  *  I-  j 

For  the  tangential  compliances  Sj^2»  ^23*  there  does  not 

seem  to  be  a  similar  rule.  bounds  do  generate  hounds,  but  It 

Is  not  known  beforehand  whether  lower/upper  C. .  bounds  will  generate 

* 

lower /upper  bounds,  or  vice  versa. 

It  la  worthwhile  to  show  how  Poisson's  ratio  bounds  are  obtained, 
Assume,  for  example,  that  and  Sj^2  hounds  are  known  by  numerical  cal¬ 
culation.  Thus, 


The  form  (15  a) 


eJ(-) 


g*(-) 

^11 


'12 


<  s*  < 

-  ^11  -  ^11 


s*(-)  <  sU  < 


_  ^12  -  ‘*12 


<  E?  < 


_*(+)  -  *'1  -  _*(-) 
^11  ^11 


-  F*W 


(183) 

(184) 

(185) 


which  yields  the  E*  bounds.  From  (16  a)  assuming  positive  Ej^, 

-  eJ  VJ2  <  -eJ  .  (186) 


If  the  left  and  right  sides  are  positive,  which  Is  mostly  the  case,  then 
because  of  (185) 

-  eJ^"^  VJ2  <  -E2^+^  S*^"^  .  (187) 

Thus,  (187)  gives  bounds  on  a  Poisson's  ratio. 

Note  that  In  a  few  cases,  negative  values  of  Poisson's  ratio 
were  obtained  by  numerical  calculation  (see  Tables  V  to  VII). 
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It  Is  of  utmost  Importance  to  compare  the  numerical  values  ob¬ 
tained  In  this  report  to  experimental  results.  Unf jrtunately,  consistent 
experimental  results  do  not  seem  to  be  available  In  the  literature. 

A  check  with  some  preliminary  experimental  results  (privately 
communicated)  was  given  In  Progress  Report  No.  8  under  this  contract. 
However,  as  pointed  out  there,  the  experimental  results  are  Inconsistent. 
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Table  I 

LOWER  BOUNDS  OF  EFFECTIVE  MODULI,  C*^  (10®  psi) 
BORON/EPOXY,  6^  =  ♦  =  30* 


VF 

C  * 

hi 

c* 

^22 

C* 

^33 

C  * 

^12 

C* 

h3 

c* 

^23 

0.1 

4.06398 

0.92263 

1.27007 

0.47054 

1.51607 

0.47850 

0.2 

7.33979 

1.07215 

1.76499 

0.51625 

2.60366 

0.52991 

0.3 

10.63586 

1.26332 

2.29864 

0.57258 

3.69748 

0.58937 

0.4 

13.96205 

1.51632 

2.88984 

0.64528 

4.79919 

0.66262 

0.5 

17.33581 

1.86696 

3.57193 

0.74467 

5.91383 

0.76029 

0.6 

20.79132 

2.38527 

4.41018 

0.89091 

7.05048 

0.90310 

0.7 

24.40579 

3.22934 

5.55201 

1.12925 

8.22932 

1.13755 

0.8 

28.39842 

4.84713 

7.41578 

1.58744 

9.50610 

1.59548 

0.9 

33.73455 

9.20441 

11.84231 

2.82342 

11.11122 

2.86432 

1.0 


63.88867  63.88846  63.88848 


15.97217 


15.97221 


15.97224 


•f 


TABLE  II  *  f. 

UPPER  BOUNDS  OF  EFFECTIVE  MODULI,  C.j  (lO^psi) 
BORON/EPOXY,  ,  ^  =  30° 


VF 

C* 

hi 

c* 

^22 

c* 

''33 

r  * 

Li2 

c  * 

''13 

r  * 

'23 

0.1 

4.06422 

0.92646 

1.27223 

0.47150 

1.51679 

0.48137 

0.2 

7.34084 

1.08908 

1.77451 

0.52049 

2.60717 

0.54260 

0.3 

10.63848 

1.30515 

2.32217 

0.58304 

3.70532 

0.62074 

0.4 

13.96716 

1.59803 

2.93580 

0.66571 

4.81451 

0.72391 

0.5 

17.34466 

2.00846 

3.65153 

0.78005 

5.94036 

0.86641 

0.6 

20.80573 

2.61589 

4.53990 

0.94856 

7.09372 

1.07606 

0.7 

24.42891 

3.59916 

5.76003 

1.22170 

8.29866 

1.41491 

0.8 

28.43675 

5.46029 

7.76067 

1.74073 

9.  07 

2.05535 

0.9 

33.80588 

10.34575 

12.48432 

3.10875 

11.32522 

3.72024 

1.0 

63.88867 

63.88840 

63.88845 

15.97215 

15.97220 

15.97219 

TABLE  III 


EFFECTIVE  MODUL;  C*.  (lO^si) 
BORON/EPOXY .  6^  =  .  ^  =  30° 


VF 

- » 

<=11 

it 

C22 

* 

^33 

* 

^12 

* 

^13 

0* 

^23 

0.1 

4.06406 

0.92390 

1.27079 

0.47086 

1.51631 

0.47946 

0.2 

7.34014 

107780 

1.76816 

0.51767 

2.60505 

0.53414 

0.3 

10.63673 

1.27726 

2.30649 

0.57607 

3.70009 

0.59983 

0.4 

13.96376 

1.54356 

2.90516 

0.65209 

4.80430 

0.68305 

0.5 

17.33876 

1.91413 

3.59846 

0.75646 

5.92268 

0.79566 

0.6 

20.79613 

2.46214 

4.45342 

0.91012 

7.06489 

0.96076 

0.7 

24.41350 

3.35262 

5.62135 

1.16006 

8.25243 

1.23000 

0.8 

28.41120 

5.05152 

7.53074 

1.63854 

9.54443 

1.74877 

0.9 

33.75833 

9,58485 

12.05632 

2.91853 

11.18255 

3.14957 

1.0 

63.88867 

63.88844 

63.88847 

15.97216 

15.97220 

15.97222 

TABLE  IV 

EFFECTIVE  MODULI,  C*.  (lO^psi) 

*  J 

BORON/EPOXY,  6^  •=  ,  ♦  =  30“ 


VF 

ir 

^11 

C22 

^33 

* 

^12 

★ 

^13 

Q  * 

^23 

0.1 

4.06414 

0.92518 

1.27151 

0.47118 

1.51655 

0.48042 

0.2 

7.34049 

1.08344 

1.77134 

0.51908 

2.60611 

0.53837 

0.3 

10.63760 

1.29120 

2.31433 

0.57956 

3.70271 

0.61028 

0.4 

13.96546 

1.57079 

2.92048 

0.65890 

4.80941 

0.70348 

0.5 

17.34171 

1.96129 

3.62499 

0.76826 

5.93152 

0.83104 

0.6 

20.80093 

2.53901 

4.49666 

0.92934 

7.07930 

1.01841 

0.7 

24.42120 

3.47589 

5.69069 

1.19088 

8.27554 

1.32245 

0.8 

28.42397 

5.25590 

7.64571 

1.68964 

9.58275 

1.90206 

0.9 

33.78210 

9.96530 

12.27032 

3.01364 

11.25389 

3.43491 

1.0  63.88867 


63.88842  63.88846  15.97216  15.97220  15.97221 


TABLE  V  g 

BOUNDS  OF  EFFECTIVE  COMPLIANCE,S*j.  ^  10 

BORON/EPOXY,  ^  =  30* 


* 

★ 

★ 

★ 

★ 

VF 

^11 

$22 

^33 

^12 

^13 

^23 

0.1 

0.44607 

1.35195 

1.66902 

0.06138 

-0.55506 

-0.58472 

0.2 

0.29177 

1.10945 

1.37568 

0.08746 

-0.45542 

-0.46774 

0.3 

0.21832 

0.90572 

1.11780 

0.07808 

-0.36923 

-0.36669 

0.4 

0.17006 

0.72741 

0.89274 

0.06247 

-0.29430 

-0.28181 

0.5 

0.13403 

0.57101 

0169685 

0.04679 

-0.22915 

-0.21161 

0.6 

0.10532 

0.43373 

0.52601 

0103269 

-0.17233 

-0.153890 

0.7 

0.081516 

0.31269 

0.37617 

0.020477 

-0.122472 

-0.106310 

0.8 

0.061062 

0.205066 

0.243670 

0.010028 

-0.078356 

-0.066742 

0.9 

0.042509 

0.108294 

0.125361 

0.001225 

-0.038927 

-0.033382 

1.0 

0.017391 

0.017391 

0.017391 

-0.003478 

-0.003478 

-0.003478 
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TABLE  VI 

BOUNDS  OF  EFFECTIVE  COMPLIANCE,  S?'.  x  10® 

*  J 

BORON/EPOXY.  ^  =  30° 


VF 

★ 

®11 

* 

®22 

* 

®33 

★ 

^12 

* 

®13 

* 

®23 

0.1 

0.446308 

1.355273 

1.670097 

0.060507 

-0.555551 

-0.582827 

0.2 

0.292335 

1.119950 

1.379173 

0.085023 

-0.456826 

-0.461687 

0.3 

0.219223 

0.924096 

1.124055 

0.074002 

-0.371605 

-0.355973 

0.4 

0.171229 

0.751987 

0.901319 

0.057120 

-0.297459 

-0.267286 

0.5 

0.135325 

0.598766 

0.706788 

0.40800 

-0.232733 

-0.194998 

0.6 

0.106599 

0.461206 

0.536071 

0.026785 

-0.175903 

-0.137265 

0.7 

0.082613 

0.336601 

0.385070 

0.015353 

-0.125597 

-0.091723 

0.8 

0.061857 

0.222670 

0.250267 

0.006287 

-0.080646 

-0.055966 

0.9 

0.042905 

0.117497 

0.128776 

-0.000686 

-0.40091 

-0.027775 

1.0 

0.017391 

0.017391 

0.017391 

-0.003478 

-0.003478 

-0.003478 
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TABLE  VII 

EFFECTIVE  COMPLIANCE  S*^  x  10® 
BORON/EPOXY,  ^  =  30“ 


VF 

CL 

S* 

33 

S* 

i  u 

S* 

^13 

S* 

^23 

0.1 

0.446231 

1.354162 

1.669738 

0.060798 

-0.555385 

-C. 583459 

0.2 

0.292144 

1.116406 

1.377995 

0.085845 

-0.456352 

-0.463731 

0.3 

0.218914 

0.917813 

1.121917 

0.075396 

-0.370792 

-0.359639 

0.4 

0.170824 

0.743450 

0.898339 

0.058978 

-0.296360 

-0.272331 

0.5 

0.134868 

0.588966 

0.703278 

0.042916 

-0.231467 

-0.200863 

0.6 

0.106144 

0.451354 

0.532464 

0.028904 

-0.174621 

-0.143226 

0.7 

0.082214 

0.327913 

0.381838 

0.017215 

-0.124462 

-0.097022 

0.8 

0.061566 

0.216218 

0.247849 

0.007658 

-0.079807 

-0.059916 

0.9 

0.042760 

0.114128 

0.127526 

0.000014 

-0.039665 

-0.029827 

1.0 

0.017391 

0.017391 

0.017391 

-0.003478 

-0.003478 

-0.003478 

1.0  0.017391 


TABLE  VIII 

EFFECTIVE  COMPLIANCE,  x  10^ 
BORON/EPOXY,  6^  =  ^  =  30® 


0.1  0.446155  1.353055  1.669379  0.061088  -0.555220  -0.584089 
0.2  0.291956  1.112906  1.376830  0.086658  -0.455883  -0.465749 
0.3  0.218613  0.911690  1.119833  0.076755  -0.369999  -0.363211 
0.4  0.170436  0.735264  0.895480  0.060760  -0.295307  -0.277168 
0.5  0.134437  0.579728  0.699970  0.04491  -0.230273  -0.206391 
0.6  0.105721  0.442220  0.529121  0.030868  -0.173433  -0.148752 
0.7  0.081850  0.319972  0.378884  0.018916  -0.123424  -0.101865 
0.8  0.061302  0.210378  0.245661  0.008899  -0.79047  -0.063491 
0.9  0.042629  0.111075  0.126393  0.000647  -0.039279  -0.031688 

0.017391  0.017391  -0.003478  -0.003478  -0.003478 


1.0  0.017391 


T 


TABLE  IX 

YOUNG'S  MODULUS, £3(10®  psi) 
BORON/EPOXY,  ^  =  30“ 


VF 

Upper  Bound 

G  =  G 

Gt  =  G,''> 

G  =  6  ^ 

Lower  Bound 

G  =  G 

0.1 

0.5990 

0.5992 

0.598R 

0.5989 

0.2 

0.7263 

0.7269 

0.7251 

0.7257 

0.3 

0.8930 

0.8946 

0.8896 

0.8913 

0.4 

1.1167 

1.1202 

1.1095 

1.1132 

0.5 

1.4286 

1.4350 

1.4149 

1.4219 

0.6 

1.8899 

1.9011 

1.8654 

1.8781 

0.7 

2.6393 

2.6583 

2.5969 

2.6189 

0.8 

4.0707 

4.1039 

3.9957 

4.0347 

0.9 

7.9118 

7.9770 

7.7654 

7.8415 

1.0 

57.5000 

57.5000 

57.5000 

57.5000 

13.  DETERMINATION  OF  BOUNDS  ON  THE  SHELL  STIFFNESS 
TENSOR  FOR  BIAXIALLY  REINFORCED  MATERIALS 


GENERAL  THEORY 


We  now  derive  expressions  for  the  components  of  the  shell 
stiffness  tensor  In  terms  of  the  moduli  for  blaxlally  reinforced  mater¬ 
ial.  It  will  be  shown  that  bounds  on  the  moduli  correspond  to  bounds 
on  the  various  components  of  the  shell  stiffness  tensor.  The  Initial 
discussion  will  be  presented  In  tensor  form,  which  Is  particularly  con¬ 
venient,  and  the  final  results  for  orthotropic  shells  of  revolution  will 
be  presented  In  terms  of  appropriate  physical  components. 


For  shell  theory,  we  use  surface  tensors,  and  components  are 
denoted  by  Greek  Indices  which  take  the  values  of  1,  2. 


The  constitutive  relations  for  shells  based  on  the  Klrchhoff 
hypothesis  take  the  form: 


-00  ^oBvA 

N  -  D  p  .  +  D,  q  .  , 

o  1  ^vX* 

oB  vXoB  „  .  ^vXoB  „ 

”  ’  °1  PvX  +  %X» 


(188) 

(189) 


^  Ol  6  01 6 

where  N  and  M  are  suitably  defined  symmetric  stress-resultants  and 
stress  couples,  respectively;  p^^  and  are  the  strain  measures  of  shell 

theory  defined  below;  and  the  quantities  (n  ■  0,1,2)  are  the  com¬ 

ponents  of  the  shell  stiffness  tensor.  The  shell  strain  measures  are 
given  In  terms  of  the  components  of  displacement  by  the  relations 


In  equations  (190)  and  (191),  are  the  membrane  displacements,  while  w 
is  the  normal  displacement;  B^  represents  the  components  of  the  rotation 
of  the  normal  to  the  middle  surface,  and  the  notation  () n  indicates 
covariant  differentiation  with  respect  to  the  surface  metric.  The  quanti¬ 
ties  b  „  and  b°  are  coefficients  of  the  second  fundamental  form  of  the 

01  p  P 

middle  surface  and  are  related  to  the  principal  radii  of  curvature  of  the 
surface. 
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By  virtue  of 
rotation  components  3^ 


the  neglect  of  transverse  shear  deformation,  the 
are  related  to  the  displacements  by  the  relations 


6 

a 


3w  .V 

T - b  V  . 

dx  a  V 
a 


(192) 


The  shell  stiffness  tensors  appearing  In  Equations  (188)  and  (189) 

are  defined  by  considering  the  strain  energy  function  W  for  the  shell  and 
requiring  that  the  following  relation  be  satisfied: 


jaB  .  A  _  h“8  .  ^  . 


(193) 


This  Insures  that  the  constitutive  relations  of  shell  theory 
are  consistent  with  the  laws  of  thermodynamics.  The  shell  strain  energy 
density  function  Is  then  written  In  the  form 

By  virtue  of  the  symmetry  of  the  stress  resultants  and  couples  and  of  the 
strain  measures,  as  well  as  the  positive  definite  character  of  W,  It  can 
be  shown  that  the  various  stiffness  tensors  satisfy  the  following 
symmetries : 


dOBY'S  .  ,  jja66Y  ^  jjY6a6 

o  o  o  o  ’ 

oOieYi  8oy6  a66Y  ,  j.Y'SaB 
^2  ^^2  *^2  ■  *^2  ’ 

d“®Y«  .  D6otY6  ,  pOBdY  ,  (I95) 

It  should  be  noted  that  the  tensor  does  not  exhibit  symmetry  with 

respect  to  the  pairs  of  Indices  (oB)  and  (y5). 

The  derivation  leading  to  the  definition  of  the  stiffness  ten¬ 
sor  can  be  found  In  [9]  and  will  not  be  repeated  here.  The  components 
of  given  In  the  equation: 
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In  Equations  (196)  through  (198),  the  various  quantities  are  defined  as 


gVXe< 

0 


h/2 

-h/2 


tvXck  .  3 
uA  dx  • 


(199) 


h/2 

(x^)”  dx^  (n  -  0,1,2).  (200) 

IT  P 

-h/2 


In  these  equations,  we  have  Introduced  the  shifters  defined  as 

OL  .a  3,  a 

“s  •  ■'b  -  ”  \  • 

while  the  "inverse  shifters"  are  defined  by  the  equation 


(201) 


and  we  have  used  the  notation 


a-e 

^6% 


(202) 


(203) 


w  "  y 


6' 


tvXck 


In  Equations  (199)  and  (200),  the  quantities  A  are  the  shifted  com¬ 
ponents  of  the  elasticity  tensor  for  the  material.  The  term  "shifted" 
Indicates  that  these  components  are  associated  with  a  system  of  base 
vectors  consisting  of  the  middle  surface  base  vectors  and  a  unit  vector 
normal  to  the  middle  surface. 


It  should  be  noted  that  the  Integrations  Indicated  In  Equations 

(199)  and  (200)  can  be  carried  out  exactly,  even  though  the  shifters 

Uq  are  functions  of  the  thickness  coordinate  x'’.  However,  for  most  cases 

^  h 

of  practical  interest.  It  Is  permissable  to  neglect  terms  of  order  — 

In  comparison  to  unity,  where  h  denotes  the  shell  thickness  and  R  the 
least  radius  of  curvature.  The  result  of  such  a  procedure  is  a  first- 
order  shell  theory  which  Is  sufficiently  accurate  for  most  cases.  In  the 
remainder  of  the  present  work,  we  will  restrict  our  consideration  to  this 
first-order  theory.  Furthermore,  to  simplify  the  results  we  assume  that 
the  layers  are  located  symmetrically  with  respect  to  the  shell  reference 

surface.  Introducing  the  approximation  -^  <<  1  In  Equations  (196)  through 

(200)  leads  to  the  following  results  for  a  single-layer  shell: 
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„»6i,s  .  joevs 

0  . 


0 

,0376 


^qi6y6  TaeY<5  ill 

“2  ~  ^  12  • 

Thus,  the  constitutive  relations  are  now  reduced  to  the  form 

'oB  r.“6Y6  toByS  . 

N  -D  p.-A  hp-. 


(204) 


M  ■  D  p  ,  “  A  TT 


y6 


12  Sfi 


(205) 


It  Is  clear  from  Equations  (204)  that.  If  we  have  bounds  on  the  components 

of  the  elasticity  tensor  we  have  corresponding  bounds  on  the  shell 

stiffness  tensor.  The  fact  that  we  are  using  the  shifted  components  of 
the  elasticity  tensor  does  not  alter  this  property,  since  the  "shifting" 
operation  depends  only  on  the  geometry  of  the  coordinate  system  chosen. 

CONSTITUITIVE  RELATIONS  FOR  ORTHOTROPIC  SHELLS 

For  an  orthotropic  material  with  planes  of  elastic  symmetry 
coinciding  with  the  coordinate  planes,  the  stress-strain  relations  may 
be  written  In  the  form 


11 

0 

^11 

^22 

^33 

0 

0 

0 

11 

e 

22 

0 

^11 

22 

33 

0 

0 

0 

22 

e 

j 

33 

0 

33 

"n 

33 

r 

^22 

33 

0 

0 

0 

33 

e 

P 

12 

0 

- 

0 

0 

0 

12 

0 

0 

12 

e 

1 

23 

0 

0 

0 

0 

c“ 

23 

0 

1 

31 

0 

0 

0 

0 

0 

0 

^31 J 

31 

e 

(206) 

the  physical 

components 

of 

the  stress. 

strain,  aid  elasticity 

ten- 

sors  have  been  used. 


Thus,  the  quantities  are  precisely  the  moduli 


previously  bounded  for  the  blaxlally  reinforced  material.  Since  the 
classical  theory  of  shells  neglects  the  effect  of  transverse  normal  stress 
and  tranverse  shear  deformation,  we  are  only  Interested  In  those  parts 
of  Equation  (206) which  refer  to  the  1  and  2  directions.  Therefore,  for  our 
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purposes,  we  may  write  a  restricted  form  of  Equation  (206)  as  follow.: 


,11 

'll 

^22 

,22 

p22 

11 

^12 

0 

0 

0 

®11 

0 

®22 
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^12_ 

_^33_ 

(207) 


We  note  that  the  moduli  In  Equations  (206)  and  (207)  satisfy  the  relations 


Now,  the  shifted  components  of  the  elasticity  tensor  In  orthog¬ 
onal  surface  coordinates  are  given  by 


70088 

„ao  88 

88 

A 

a 

ja6a8 

-  c““  a““ 
08 

a^® 

(no  sum  on  a, 8)  > 
(no  sum  on  a, 8)  , 


(208) 


where  a^^  are  the  components  of  the  surface  metric  and  the  a 


a8 


are  the 


corresponding  components  of  the  conjugate  tensor,  defined  such  that 


a  o  ® 
a8 


-  .7 


(209) 


We  note  that  for  orthogonal  surface  coordinates,  we  have  a^^  -0  (of*  8), 


and  It  follows  that 

a®®  ■  — ^  (no  sum).  (210) 

a 

aa 

In  order  to  Identify  the  various  components  of  the  shell  stiffness  tensor 
In  terms  of  the  previously  bounded  moduli  for  blaxlally  reinforced  mater¬ 
ial,  It  Is  convenient  to  rewrite  Equation  (205)  In  matrix  form  In  terms 
of  the  physical  components  of  the  various  tensors.  When  this  Is  done, 
we  note  that  for  an  orthotropic  material  where  planes  of  symmetry  coincide 
with  the  lines  of  curvature  of  the  middle  surface,  the  stress-strain 
relations  may  be  separated  Into  two  groups  Involving  normal  and  shear 
components,  respectively.  The  resulting  equation  may  then  be  written  In 
the  form 
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where  the  sjnnbol  <>  denotes  physical  components. 

Substituting  Equation  (204)  into  Equations  (211)  and  (212)  and 
making  use  of  Equations  (208)  through  (210)  yields  the  required  stress- 
strain  relations  In  terms  of  the  previously  bounded  moduli;  the  result  Is 
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We  see  from  Equations  (213)  and  (214)  that  we  need  to  evaluate  the  fol¬ 
lowing  stiffnesses: 

„1L  ^11  h^  ^22.  -11. 

^22  p22,  _22  hi  -12  „12  . 

^11  12  “  ^22  12  »  ^22  12*  ^12^»  ^12  12 


In  terns  of  the  slx-by-slx  matrix  notation  employed  In  previous  chapters, 

n  A 

the  moduli  c“  are  Identified  from  the  equation 
y4 
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(216) 

It  Is  now  clear  from  Equation  (216)  that  bounds  on  the  shell  stiffness 

tensor  are  obtained  from  bounds  on  the  moduli  (1,  J  -  1...6)  as 

previously  obtained  In  Section  11.  We  note  from  equations  (213)  through 
.  *  *  *  * 

(216)  that  If  bounds  on  the  quantities  C^2»  ^22*  ^44  known, 

corresponding  bounds  on  all  nonvanishing  components  of  the  shell  stiff¬ 
ness  tensor  are  Immediately  known. 


14.  ALTERNATE  SOLUTION  FOR  BOUNDS  ON 

COMPONENTS  OF  THE  SHELL  STIFFNESS  TENSOR 


In  Section  13.  bounds  on  the  components  of  the  shell  stiffness 
tensor  were  obtained  In  terms  of  the  bounds  previously  determined  for  the 
moduli  of  blaxlally  reinforced  materials. 


This  result  was  obtained  by  considering  the  shell  as  homogeneous 
and  constructed  from  a  material  having  the  appropriate  moduli  for  a  bl¬ 
axlally  reinforced  material.  In  the  present  section,  we  present  an  al¬ 
ternate  approach  based  on  the  use  of  a  layered  shell  model  hereafter 
designated  by  LSM. 


In  the  LSM  analysis,  the  shell  Is  considered  as  constructed 
from  m  layers  of  unlaxlally  reinforced  material  with  fibers  oriented  at 
angles  (p^  (1  ■  1,2, - m)  with  respect  to  the  lines  of  curvature  of  the 

shell  reference  surface.  We  again  consider  a  classical  shell  theory 
under  the  Klrchhoff  hypothesis  and  limit  the  analysis  to  the  first-order 

shell  theory  (—  <<  1).  Under  these  conditions,  the  components  of  the 

shell  s*'lffness  tensors  are  given  by 


DaSyS 
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(n  -  0,2)  , 


(217) 


where  ^  Is  the  distance  of  any  point  from  the  reference  surface,  de¬ 
notes  the  distance  from  the  reference  surface  to  the  Inner  surface  of  the 
til  “  oiSy^ 

1  layer,  and  A.  ^  are  shifted  components  of  the  elasticity  tensor  for 
th  ^ 

the  1  layer.  When  Equation  (217)  Is  transformed  to  physical  components, 
we  obtain  the  following  results 
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(219) 


In  Equations  (218)  and  (219),  the  quantities  C.  are  the  elastic  moduli 

of  the  layer  referenced  to  the  Xj^,  X2,  coordinate  system  of  Figure  1. 

These  moduli  are  those  of  a  unlaxlally  reinforced  material  after  trans¬ 
formation  to  the  proper  coordinate  system  and  are  Identical  to  the 
quantities  or  Equations  (75)  and  (76), 
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It  Is  convenient  to  express  Equations  (218)  and  (219)  In  tesms 
of  the  thickness  h^  of  the  1^^  layer  and  the  coordinate  which  will 
represent  the  distance  from  the  reference  surface  to  the  middle  surface 
of  the  1^^  layer.  In  terms  or  these  quantities,  (218)  and  (219)  become 
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Equations  (220)  and  (221)  may  be  used  to  calculate  the  shell  stiffness 
tensor  for  the  LSM  for  any  configuration  of  layers,  provided  that  the 
transformation  of  coordinates  for  each  layer  Is  performed  using  the 
appropriate  angle  (|)^  In  the  transformation  equations  of  Section  6.  We 

note  that  the  shell  will  exhibit  complete  anisotropy  In  the  general  case. 


We  now  consider  the  special  case  of  a  blaxlally  reinforced  shell 
having  an  even  number  of  layers.  We  assume  that  the  layers  are  arranged 
on  either  side  of  the  shell  middle  surface  so  that  If  the  layer  with 
middle  surface  at  has  a  fiber  angle  4)^,  the  layer  with  middle  surface 

at  has  a  fiber  angle  (180-<|).),  IJnder  these  conditions, the  elastic 

^  ^  « 0  ft 

moduli  for  the  two  layers  become  and  "C.  ,  which  satisfy  the  sym- 

metry  conditions  (78) .  Because  of  symmetry  of  Equations  (220,  221)  with 
respect  to  we  can  then  express  the  shell  stiffness  tensors  as 
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Note  that  Equation  (224)  Is  Identical  to  (79).  We  note  that,  because 
of  the  symmetries  outlined  In  Section  6,  the  only  nonvanishing  components 


of  5,**^  are 

S6 

-  11  -  22  -  11  _  g  22  -  12 

1  11*  S  22*  ^1  22  ^1  11*  12 


(225) 


Recall  that  the  Greek  Indices  do  not  take  on  the  value  of  3.  Therefore, 
this  type  of  construction  gives  rise  to  an  orthotropic  shell,  as  should 
be  expected  from  the  results  of  Sections  11  and  13.  The  constitutive 
relations  for  shells  based  on  the  LSM  are  then  written  In  a  form  ana¬ 
logous  to  (213)  and  (214) : 


(226) 


(227) 


where  the  stiffness  components  are  given  In  Equations  (222)  through  (224). 
In  the  notation  of  Section  6  (see  C31)),  the  moduli  appearing  In  (225)  can 
be  Identified  In  the  six-by-six  matrix  notation  according  to  scheme 

C  -  C  11,  2  22  ■  S  22,  ^  22  "  ^  12*  ^  12  "  ^44  »  ^228) 

and  the  quantities  Cij  are  given  by  Equations  (82)  In  terms  of  the  single¬ 
layer  moduli  of  Equations  (42)  and  (43).  Since  four  of  these  moduli 

are  known  exactly  and  the  fifth  Is  bounded  as  a  result  of  the  work  of 
Hashln  and  Rosen  [1],  It  Is  clear  that  an  argument  similar  to  that  of 
Section  11,  leading  to  (175)  can  be  applied  to  obtain  bounds  for  the 
shell  stiffness  tensors  as  defined  by  the  LSM. 

Therefore,  we  have  succeeded  In  obtaining  two  different  sets  of 
bounds  for  the  shell  stiffness  tensors.  It  Is  likely  that  both  of  these 
sets  of  bounds  are  In  less  agreement  because  of  the  manner  In  which  the 
moduli  bounds  for  the  blaxlally  reinforced  material  were  developed. 
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We  note  that  the  results  obtained  in  Section  13  and  14  proceed 
from  the  use  of  a  first-order  shell  theory  based  on  the  Kirchhoff  hypo¬ 
thesis.  Therefore,  the  transverse  shear  deformation  and  transverse  nor¬ 
mal  stress  have  been  neglected,  and  the  application  should  be  limited  to 
sufficiently  thin  shells.  When  applied  to  such  shells,  these  results 
should  yield  acceptable  results  for  deformations  as  well  as  overall 
stress  resultants  and  stress  couples.  Certain  quantities  such  as  the 
interlayer  shear  stresses  cannot  be  calculated  directly  from  such  a 
theory  but  can  be  determined  approximately  after  the  displacement  field 
is  known.  Thus,  the  results  presented  in  Sections  13  and  14  are  subject 
to  the  approximations  of  a  first-order  classical  shell  theory  and  are 
consistent  with  these  approximations. 
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15.  CONCLUSION 


Bounds  for  the  orthotropic  effective  elastic  ccdull  and 
compliances  of  blaxlally  fiber  reinforced  materials  have  been  derived  by 
variational  procedures.  It  has  been  assumed  that  the  single-layer  unl- 
axlally  reinforced  layers  are  described  by  the  composite-cylinder  assem¬ 
blage  model,  [1]. 

Numerical  evaluation  of  the  bounds  showed  that  they  are  ex¬ 
tremely  close  together  In  the  majority  of  cases.  Thus,  they  should  pro¬ 
vide  good  estimates  for  elastic  properties  of  blaxlally  reinforced 
materials.  The  results  obtained  for  effective  elastic  properties  have 
also  been  exploited  to  construct  bounds  for  the  membrane  and  bending 
stiffnesses  of  orthotropic  blaxlally  reinforced  shells. 

It  Is  of  utmost  Importance  to  compare  the  present  theoretical 
results  with  experimental  findings.  Unfortunately,  consistent  experi¬ 
mental  results  for  elastic  properties  of  blaxlally  reinforced  materials 
do  not  seem  to  be  available  In  the  literature.  There  Is  accordingly 
great  need  for  such  an  experimental  program. 
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